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surfaces,  one  a  reference  plane  and  the  other  an  effective  plane  created  by  the 
dielectric  d;scontinuity  at  the  burning  propellant  surface.  In  the  pressure  range  from 
500  to  10,000  psig,  steady-state  burning  rate  tests  with  a  single  base  nitrocellulose  gun 
propellant,  M6,  were  successfully  conducted  in  a  large  volume  combustion  fixture 
which  incorporated  a  coaxial-to-waveguide  adapter  having  unique  pressure  sealing 
microwave  feedthroughs.  In  addition  to  yielding  burning  rate  values  intermediate  to 
those  from  conventional  closed  chamber  and  strand  burner  tests,  the  experiments 
established  the  uncertainty  in  the  proportionality  constant  linking  M6  burning  rate  and 
the  time  derivative  of  relative  phase  angle.  A  one-dimensional  theoretical  model  was 
formulated  and,  with  the  aid  of  appropriately  chosen  constants,  yielded  predictions 
which  agreed  very  well  with  the  experimentally  observed  behavior  of  a  solid  propellant 
undergoing  planar  burning  within  a  waveguide..  In  conjunction  with  a  nonlinear 
parameter  estimation  technique,  utilization  of  the  model  allowed  apparent  instantaneous 
burning  rates  to  be  corrected  for  biases  which  had  been  caused  by  the  presence  of 
significant  fixed  microwave  reflections. 

In  the  attempt  to  extend  the  microwave  technique  to  rapid  pressurizations,  a  twin 
propellant  strand  differential  system  was  constructed  to  help  isolate  the  propellant 
combustion  response  from  its  mechanical  response.  In  this  approach,  changes  in  relative 
phase  angle  were  intended  to  be  associated  with  only  that  motion  of  the  microwave 
reflecting  plane  caused  by  propellant  burning.  The  ability  of  the  technique  to  cancel 
propellant  mechanical  responses  was  evaluated  in  a  series  of  vibration  or  non-burning 
dynamic  tests  conducted  at  two  distinctly  different  pressurization  rates,  slow  and  fast. 
Within  the  volume  of  interest,  typical  pressure  changes  were  from  1450  to  7000  psig 
and  were  created  by  a  controlled  influx  of  nitrogen  gas;  characteristic  pressurization 
rates  were  50  X  103  and  1  X  106  psi/sec.  Dynamic  experiments  which  included 
propellant  combustion  were  also  conducted  under  similar  pressurization  conditions. 
Results  from  the  transient  burning  tests  and  the  vibration  experiments  indicated  that 
major  problems  must  still  be  overcome  before  accurate,  direct  dynamic  burning  rate 
measurments  can  be  made  via  the  microwave/Doppler  phase  shift  technique. 
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NOMENCLATURE 


A  =  magnitude  of  a  vector  representing  a  microwave  fixed  reflection 

A  =  cross  sectional  area  of  one  rectangular  opening  in  the  dual 

propellant  waveguide 

Am  =  relative  amplitude  based  on  the  DFT  amplitude  of  a  given 

physical  quantity;  used  in  spectral  analysis  of  experimental  data 

Ap  =  cross  sectional  area  of  the  solid  propellant  strand 

Are,  =  change  in  decibels  between  the  test  and  reference  RF  voltage 

amplitudes  on  a  dynamic  test 

A,  =  throat  area  in  a  flow  nozzle 

a  =  exponent  for  pressure  as  used  in  a  steady  burning  rate  correlation 

a  =  dimension  of  the  longest  side  in  the  cross  section  of  a 

rectangular  waveguide;  used  in  defining  cutoff  wavelength 

a  =  length  of  X-band  waveguide  as  per  Figure  9 

B  =  magnitude  of  a  vector  representing  the  microwave  reflection  from 

a  moving  surface  which  is  perpendicular  to  the  incident  wave  train 

B  =  coefficient  in  a  steady  burning  rate  correlation 


xvii 

Bu  =  constant  associated  with  the  initial  magnitude  of  vector  B 

b  =  length  of  X-band  waveguide  as  per  Figure  9 

bj|  =  discrete  filter  weights  for  differentiating 

b*  =  discrete  filter  weights  for  smoothing 

C  =  magnitude  of  a  vector  resulting  from  the  sum  of  vector  A  and 

— ► 

vector  B;  on  steady  burning  tests  C  is  the  test  signal  magnitude  as 
measured  by  the  microwave  technique 

Cd  =  discharge  coefficient 

C  .  =  any  one  of  a  set  of  discrete  measured  test  signal  magnitudes 

C,  j  =  any  one  of  a  set  of  theoretically  predicted  test  signal  amplitudes 

c  =  length  of  X-band  waveguide  as  per  Figure  9 

c  =  speed  of  light 

cp  =  specific  heat  at  constant  pressure 

cy  =  specific  heat  at  constant  volume 

D  =  flexural  rigidity 

D  =  diameter 

d,  =  flow  nozzle  throat  diameter 
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=  total  effective  energy  available  from  a  unit  mass  of  propellant 

=  elastic  modulus 

=  accumulative  energy  based  on  the  DFT  energy  of  a  given 
experimental  quantity;  used  in  spectral  analysis  of  experimental  data 

=  specified  tolerance  or  relative  error 

=  internal  energy  per  unit  mass 

=  independent  values  of  the  discrete  Fourier  transform  or  DFT 
=  frequency 

=  cutoff  frequency  for  the  TE,0  wave  in  a  dielectric  filled 
rectangular  waveguide 

=  any  frequency  within  the  range  of  nominal  system  bandwidths 
=  reference  frequency  in  the  range  of  nominal  system  bandwidths 
=  sampling  frequency  of  the  A/D  system 

=  shear  force 

=  end  of  the  pass  band  for  digital  filtering 
=  beginning  of  the  rejection  band  for  digital  filtering 


=  shear  modulus 
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H(z) 


H(oj) 


H(co) 


h 

h 

kf 
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k 

ki 
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=  Newton’s  proportionality  constant 

=  individual  members  of  the  differentiated  output  set 

=  individual  members  of  the  smoothed  output  set 

=  Heaviside  unit  function;  equal  to  1  for  positive  argument  and  0 
for  negative  argument 

=  real  symmetric  transfer  function  associated  with  an  infinite 
number  of  filter  weights 

=  real  symmetric  transfer  function  associated  with  a  finite  number 
of  filter  weights 

=  plate  thickness 

=  enthalpy 

=  force  or  force  constant  of  propellant 

=  wave  number 

=  constant  used  in  defining  the  size  of  the  transition  band  for  the 
digital  filter 

=  effective  spring  constant 
=  spring  constant  based  on  EA/h 

=  spring  constant  based  on  the  deflection  of  a  circular  flat  plate 
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£  =  instantaneous  length  of  a  propellant  strand  affixed  within  a 

rectangular  waveguide 

6  =  dft/dt;  instantaneous  velocity  of  the  microwave  reflecting  plane 

associated  with  the  exposed  surface  of  a  solid  propellant  strand 
affixed  within  a  rectangular  waveguide 

2o  =  initial  length  of  propellant  strand 

m  =  mass 

mp  =  efflux  of  propellant  combustion  products  in  mass  per  unit  time 

N  =  number  of  discrete  samples  of  data 

N  =  associated  with  the  number  of  filter  weights;  total  number  of 

filter  weights  equals  2N  +  1 

n  =  number  of  moles 

P  =  power 

p  =  pressure 

q  =  uniformly  distributed  force  per  unit  area 

R  =  universal  gas  constant 

Ro  =  radius  of  a  circular  flat  plate 

R  =  gas  constant  for  a  particular  gas 
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r  =  propellant  burning  rate 

r  =  pressure  ratio;  outlet  or  low  divided  by  inlet  or  high  static 

pressure 

r  =  radius 

rcrit  =  the  ratio  between  the  outlet  and  inlet  static  pressures  at  throat 
sonic  velocity 

T  =  time  interval  between  samples 

T  =  temperature 

To  =  isochoric  adiabatic  flame  temperature 

Tref  =  reference  temperature 

t  =  time 

t  =  transformation  variable 

V  =  volume 

V  =  voltage 

V  =  magnitude  of  the  vector  representing  the  microwave  reflection 
from  the  reference  plane  on  a  steady-state  burning  test 

vm  =  mechanical  response  velocity  of  the  surface  of  the  reference 

propellant  strand 
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vm  t  =  mechanical  response  velocity  of  the  surface  of  the  test  propellant 
strand 

W  =  mass  rate  of  flow  through  a  nozzle 

w  =  local  deflection  of  a  circular  flat  plate 

w  =  mean  deflection  of  a  circular  flat  plate 

X  =  dial  setting  of  the  digital  test  channel  gain  control  on  the 

network  analyzer 

Xnul|  =  test  channel  gain  in  dB  which  if  internally  imposed  by  the 

network  analyzer  would  make  the  output  from  the  relative  gain 
indicator  equal  to  zero  volts 

x  =  mass  fraction 

xbfW  =  individual  values  of  the  components  in  the  new  base  vector 

|xnew|.  each  component  relates  to  an  Xj  coordinate  axis 

x°L°  =  individual  values  of  the  components  in  the  old  base  vector 

Jxold|.  eac|j  component  relates  to  an  x(  coordinate  axis 

xp  j  =  individual  values  of  the  components  in  the  vector  |xpj  which 
evolves  from  a  pattern  move;  each  component  relates  to  an  x; 
coordinate  axis 

Xj  =  coordinate  axis;  there  are  n  coordinate  axes  or  dimensions 

Y  =  instantaneous  relative  amplitude  based  on  output  from  the  relative 

gain  indicator 
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=  expansion  factor 

=  initial  relative  amplitude  based  on  output  from  the  relative  gain 
indicator 

=  input  impedance  of  the  harmonic  frequency  converter 

=  angle  between  a  vector  representing  a  microwave  fixed  reflection 
and  the  real  axis  in  a  phasor  plot 

=  attenuation  constant  for  the  one-way  loss  of  a  TE,0  wave  in  a 
dielectric  filled  rectangular  waveguide 

=  phase  constant  for  a  propellant  filled  rectangular  waveguide 
=  phase  constant  for  standard  X-band  waveguide 
=  ratio  of  specific  heats 

=  ratio  of  specific  heats  at  the  isochoric  adiabatic  flame  temperature 

=  instantaneous  change  in  length  of  the  test  propellant  strand 

=  instantaneous  change  in  the  relative  amplitude  based  on  output 
from  the  relative  gain  indicator 

=  total  change  in  actual  phase  angle  during  complete  combustion  of 
the  test  propellant  strand 
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A\pm  =  total  change  in  measured  phase  angle  during  complete  combustion 
of  the  test  propellant  strand 

5  =  loss  angle 

5|  =  fixed  exploratory  increment  parallel  to  the  x;  coordinate  axis 

e  =  permittivity 

e  =  real  permittivity 

e"  =  complex  permittivity 

eo  =  permittivity  of  free  space 

ef  =  relative  permittivity,  commonly  called  the  dielectric  constant 

17  =  multivariable  objective  function 

T)a  =  covolume  in  Noble-Abel  equation  of  state 

X  =  nondimensional  radian  frequency 

X.  =  nondimensional  end  of  the  radian  frequency  pass  band  in  the 

digital  filter 

Xc  =  cutoff  wavelength  for  the  TE10  wave  in  a  dielectric  filled 

rectangular  waveguide 

X  =  phase  wavelength  of  the  TK)0  wave  in  a  dielectric  filled 

rectangular  waveguide 
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XR  =  nondimensional  sharpness  of  the  radian  frequency  roll-off  after 

the  end  of  the  digital  filter  pass  band 

XR  =  nondimensional  parameter  which  specifies  the  size  of  the  radian 

frequency  region  within  which  errors  between  H(co)  and  H(co)  are 
of  a  specific  magnitude 

p  =  attenuation  constant  based  on  actual  phase  angle  and  relevant  to 

the  two-way  loss  of  a  TE)0  wave  in  a  dielectric  filled  waveguide 

p'  =  attenuation  constant  based  on  propellant  test  strand  length  and 

relevant  to  the  two-way  loss  of  a  TE,0  wave  in  a  dielectric  filled 
waveguide 

p0  =  permeability  of  free  space 

v  =  Poisson’s  ratio 

v  =  frequency 

p  =  density 

pp  =  density  of  solid  propellant 

a  =  standard  deviation 

t  =  shear  stress 


r 


=  time 
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=  phase  noise  at  any  frequency  within  the  range  of  nominal  system 
bandwidths 

0ref  =  reference  phase  noise  value  in  the  range  of  nominal  system 
bandwidths 

=  actual  phase  angle;  defined  in  Figures  10  and  11 

i//a  =  di/^/dt;  time  rate  of  change  of  actual  phase  angle 

i//a  f  =  value  of  actual  phase  angle  at  propellant  burnout 

\pa  .  =  value  of  actual  phase  angle  at  the  start  of  propellant  burning 

R  =  angle  between  the  vector  BR  and  the  real  axis  in  the  phasor  plot 
of  Figure  10 

— y 

4>.d  T  =  angle  between  the  vector  BT  and  the  real  axis  in  the  phasor  plot 
of  Figure  10 

=  measured  phase  angle;  defined  in  Figures  10  and  11 
i^m  =  d\pm  I dt;  time  rate  of  change  of  measured  phase  angle 

t^m  .  =  value  of  measured  phase  angle  at  t  =  0  on  a  vibration  test 

— y 

R  =  angle  between  the  vector  CR  and  the  real  axis  in  the  phasor  plot 
of  Figure  10 

T  =  angle  between  the  vector  CT  and  the  real  axis  in  the  phasor  plot 
of  Figure  10 
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i|/SR  =  phase  shift  caused  by  a  jump  in  wave  impedance  at  the 
non-burning  propellant  surface 

\p5T  =  phase  shift  caused  by  a  jump  in  wave  impedance  at  the  burning 
propellant  surface 

t//  m  =  phase  angle  constructed  from  the  sum  of  measured  phase  angle 

and  the  change  in  theoretical  phase  angle;  relevant  only  to  vibration 
tests 

=  theoretical  phase  angle  defined  in  equation  (35);  relevant  only  to 
vibration  tests 

u>  =  radian  frequency 

cjd  =  Doppler  angular  frequency 

cjs  =  radian  sampling  frequency  of  the  A/D  system 

cjj  =  radian  frequency  at  the  end  of  the  pass  band  for  the  digital 

filter 

co2  =  radian  frequency  at  the  start  of  the  rejection  band  for  the  digital 
filter 


Subscripts 

a  =  actual 


FL 


=  flow  onset  from  a  chamber 
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fa  =  filter-amplifier;  3  dB  amplitude  point 

h  =  high  pressure  side  of  nozzle 

IG  =  ignition  of  propellant 

i  =  index  which  takes  on  the  values  1  and  2  in  order  to  identify 

parameters  relevant  to  the  driver  and  test  chambers,  respectively 

C  =  low  pressure  side  of  nozzle 

m  =  measured 

na  =  network  analyzer;  3  dB  amplitude  point 

pm  =  phasemeter;  3  dB  amplitude  point 

R  =  reference 

sys  =  nominal  system;  3  dB  amplitude  point 

T  =  test 

1  =  driver  chamber 

2  =  test  chamber 

Double  Subscripts 

i,N2  =  nitrogen  in  the  driver  and  test  chambers  for  i  =  1  and  2, 

respectively 

t. 

♦ 
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i,p  =  propellant  combustion  products  in  the  driver  and  test  chambers 

for  i  =  1  and  2,  respectively 

R,dBm  =  reference,  decibels  referred  to  1  mW 

T,dB  =  test,  decibels  referred  to  1  mW 

m 
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1.0  INTRODUCTION 


The  immediate  objectives  of  this  research  effort  were  twofold:  first,  to 
measure  the  burning  velocity  of  a  solid  propellant  responding  within  a  steady  high 
pressure  environment  via  a  microwave  Doppler/phase  shift  technique  similar  to  that 
employed  by  other  researched 1  '5 *  but  in  relatively  low  pressure  investigations;  and 
secondly,  to  further  modify  this  microwave  approach  and  evaluate  the  ability  of  the 
new  system  to  distinguish  the  combustion  response  from  the  combined 
combustion/mechanicai  response  of  a  solid  propellant  surface  when  burning  within  a 
highly  transient  pressure  environment.  The  primary  reason  for  implementing  a 
microwave  technique  was  the  potential  for  acquiring  the  one-dimensional  motion  of 
a  steady  or  dynamically  burning  propellant  strand  via  essentially  continuous,  direct 
measurement. 

In  the  case  of  transient  burning  rate  determination,  the  current  primary 
tool  is  the  closed  bomb  or  closed  chamber  test'6, 7  *  which  is  inherently  indirect 
since  it  relates  the  experimentally  observed  pressure-time  data  to  the  burning  rate  by 
an  interactive  framework  of  assumed  theoretical  expressions  for  the  thermodynamic 
and  thermochemical  processes  taking  place  in  a  closed  volume.  In  a  recent 
publication/8*  a  comparison  was  made  of  the  theory,  data  reduction  techniques,  and 
overall  methodology  utilized  by  researchers  active  in  closed  bomb  work.  The  beauty 
of  a  direct  technique  for  determining  solid  propellant  burning  rates  in  rapidly 
varying  pressure  environments  is  that  the  acquired  data  would  show  the  deviation 
from  steady-state  predictions  of  burning  behavior  and  could  be  used  to  evaluate  and 
improve  theoretical  models  of  the  instantaneous  regression  rate  of  solid  propellant 
under  rapid  pressure  transient.  Such  models  of  dynamic  burning  have  been 
numerically  compared'9'1 1  *  for  a  given  imposed  pressure  variation  at  the  edge  of 


*Parenthetical  references  placed  superior  to  the  line  of  text  refer  to  the 
bibliography. 
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the  flame  zone.  Although  strong  differences  in  the  dynamic  burning  rate  predictions 
were  observed,  no  attempt  could  be  made  to  identify  which  was  even  qualitatively 
correct  because  of  a  lack  of  existing  experimental  data. 

An  excellent  review  article  on  the  theory  of  dynamic  burning  of  solid 
propellants  has  recently  been  published. ( 1 The  dynamic  burning  effect  is  simply 
the  departure  of  the  instantaneous  burning  rate  from  the  steady-state  value 
corresponding  to  the  instantaneous  pressure.  This  effect  is  physically  caused  by  the 
finite  time  interval  required  for  the  temperature  profiles  inside  the  condensed  phase, 
and  possibly  the  reaction  zone,  to  adapt  to  the  transient  gas  phase  pressure.  When 
the  characteristic  time  of  the  pressure  transient  is  of  the  same  order  or  less  than 
the  characteristic  time  of  the  unburned  solid  phase,  dynamic  burning  behavior 

becomes  important. 

A  variety  of  direct  measurement  techniques  have  been  utilized  for  the 

acquisition  of  solid  propellant  burning  rates.  Some  of  these  approaches  were 

reviewed  in  detail  in  reference  2  and  are  briefly  noted  in  the  following  list. 

1.  Probes  are  embedded  in  the  propellant  at  given  distances  and  detect 
the  arrival  time  of  the  burning  surface.  The  sensing  elements  can  be 
thermocouples,  fuse  wires,  conducting  bridges,  or  anything  which 
detects  the  passage  of  the  burning  plane. 

2.  A  positioning  servomechanism  feeds  a  propellant  strand  into  a 

combustion  chamber  at  the  same  rate  at  which  it  is  consumed.  Various 
sensing  systems  can  be  used  to  locate  the  burning  surface  and  supply 
an  accurate  feedback  signal  for  the  servomotor. 


3.  After  leaving  a  source  and  passing  through  a  test  rocket  motor  in 
which  a  solid  propellant  grain  is  burning,  a  beam  of  x-rays  impinge 
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upon  a  fluorescent  screen  of  an  image  intensifier.  The  fluoroscopic 
x-ray  images  are  then  photographed  by  a  movie  camera. 

4.  A  propellant  strand  is  used  as  the  dielectric  material  in  a  capacitor 
which  is  part  of  the  total  capacitance  of  a  tuned  L-C  circuit  in  an 
electronic  oscillator.  The  variation  in  capacitance  caused  by  burning  is 
proportional  to  the  change  in  strand  length. 

5.  In  the  ultrasonic  pulse-echo  technique,  a  high  frequency  sound  pulse  is 
transmitted  through  an  acoustic  spacer  and  a  mating  propellant  sample, 
reflects  from  the  fixed  and  moving  interfaces,  and  returns  to  the 
original  source  where  time  dependent  electrical  signals  are  produced 
and  displayed  on  a  CRT.  High  speed  photography  of  the  screen 
coupled  with  a  linear  relationship  between  the  distance  traveled  by  the 
two  displayed  signals  allows  the  determination  of  the  actual  propellant 
length  burned. 

6.  In  the  microwave  interferometric  or  CW  radar  approach,  a  microwave 
signal  at  a  fixed  source  frequency  is  passed  through  a  propellant  strand 
and  reflects  from  the  burning  surface.  The  Doppler  shifted  reflected 
signal  is  mixed  with  a  portion  of  the  original  signal  to  produce  a 
traveling  wave  with  amplitude  maxima  which  occur  at  a  frequency 
proportional  to  burning  rate. 

All  of  the  above  mentioned  direct  measurement  techniques  give  satisfactory 
burning  rate  results  for  steady  and  near  steady-state  conditions.  In  fact,  the  classical 
method  of  determining  steady-state  propellant  burning  rates  is  approach  number  one, 
better  known  as  the  strand  burner  technique.  However,  none  of  these  approaches 
has  the  extremely  small  propellant  length  resolution  and  rapid  response  time  required 
to  obtain  regression  rate  data  under  highly  transient  pressure  conditions.  To  gain 
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some  appreciation  of  the  order  of  magnitude  of  such  stringent  length  resolution  and 
response  time  requirements,  consider  the  burning  response  of  a  typical  gun 
propellant  during  a  pressurization  from  1500  psig  to  7500  psig  over  the  course  of 
6  msec.  The  average  pressurization  rate  is  1X1 06  psi/sec  and  the  propellant  will 
regress  0.0055  in.  if  it  burns  according  to  a  steady-state  correlation  and  the 
instantaneous  pressure.  If  knowledge  of  the  strand  length  at  14  spatial  locations 
would  satisfactorily  define  the  combustion  response  during  the  pressurization,  then 
the  required  length  resolution  must  be  3.93  X  10" 4  in.  (=  10  /am)  or  less.  The 
strand  length  data  must  be  acquired  at  a  rate  of  2.33  kHz  or  higher.  Only  the 
microwave  Doppler/phase  shift  technique  exhibits  the  apparent  potential  for  meeting 
the  rigorously  exacting  response  time  and  length  resolution  requirements  fundamental 
to  direct  burning  rate  determination  during  a  rapid  pressure  transient. 
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2.0  DESIGN  RATIONALE 


To  meet  the  objective  of  acquiring  high  pressure  steady-state  burning  rates 
of  solid  propellant  via  the  microwave  technique,  the  primary  obstacle  to  overcome 
was  the  requirement  of  containing  high  pressure  gas  inside  the  combustion  fixture 
while  simultaneously  introducing  the  guided  microwaves  into  the  non-burning  base 
end  of  the  propellant  strand.  In  the  case  of  transient  burning  rate  determination, 
however,  the  problems  transcend  those  associated  with  sealing  a  large  volume 
pressure  vessel.  The  motion  of  a  burning  propellant  surface  subjected  to  a  pressure 
transient  will  be  the  result  of  both  a  combustion  response  and  a  purely  mechanical 
response.  The  magnitude  of  this  mechanical  response  depends  upon  the  geometry 
and  rheological  properties  of  the  propellant  strand,  the  sample  confinement 

characteristics,  and  the  pressurization  or  depressurization  rate.  Hence,  the  primary 
obstacle  to  progress  in  transient  burning  rate  measurements  was  initially  identified  as 
the  inability  of  the  previously  utilized  microwave  techniques  to  observe  only  the 
combustion  response. 

In  the  approach  taken  to  rectify  this  deficiency,  two  guided  microwave 

signals  oscillating  at  the  same  source  frequency  are  simultaneously  passed  through 
the  base  ends  of  two  identical  propellant  strands  and  allowed  to  reflect  from  the 
opposite  end  of  each  strand.  In  one  strand  the  signal  is  reflected  from  a  burning 
surface  and  becomes  the  test  signal.  In  the  other  strand  the  signal  is  reflected  from 
the  non-burning  surface  and  becomes  the  reference  signal.  If  the  reflecting  surfaces 
of  each  propellant  strand  are  exposed  to  the  same  pressure  environment  and  if  each 
strand  has  identical  confinement  geometry,  then  the  reflecting  surfaces  should  ideally 
undergo  the  same  mechanical  response  to  the  transient  pressure  force.  It  is,  however, 
only  the  burning  surface  which  should  experience  the  additional  motion  caused  by 
the  combustion  response.  By  continuously  comparing  the  phase  angle  of  the  test 

signal  with  the  phase  angle  of  the  reference  signal,  a  relative  phase  angle  is 
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obtained.  The  rate  of  change  of  this  relative  phase  angle  is  directly  proportional  to 
the  burning  rate  of  the  propellant. 

The  pressure  transient  affecting  the  microwave  reflecting  surfaces  of  the 
propellant  samples  was  to  be  predominantly  the  result  of  a  precise  external  mass 
addition  of  inert  gas  rather  than  a  pressure  rise  caused  by  the  deflagration  of  a 
propellant  strand.  This  approach  insures  that  the  complete  pressure  history  is  the 
result  of  an  imposed,  experimentally  controllable  action.  Consequently,  the  test 
vehicle  was  envisaged  as  a  combustion  fixture  having  two  independently  pressurizable 
chambers  connected  by  a  controllable  flow  path.  The  larger  chamber  was  identified 
as  the  driver  chamber  while  the  microwave  reflecting  propellant  surfaces  were  to  be 
within  a  small  volume  test  chamber. 

To  assist  in  establishing  the  final  geometric  configuration  of  the 
dual-chambered  combustion  fixture,  a  gas  dynamic  model  was  formulated  and 
utilized  to  predict  combustion  fixture  performance.  This  theoretical  model  is 
described  in  Appendix  C.  Use  of  the  model  allowed  parametric  investigations  of  all 
important  system  parameters;  for  example,  consideration  was  given  to  initial 
temperatures  and  pressures  and  gaseous  compositions  in  the  two  chambers,  the  size 
of  the  connecting  flow  nozzle,  chamber  initial  volumes,  propellant  sample  burning 
rate  and  cross  sectional  area,  and  ignition  and  flow  onset  times. 

One  of  the  more  important  outputs  from  the  gas  dynamic  model  was  the 
pressure  history  in  the  test  chamber.  As  shown  in  Figures  1  and  2,  there  was 
reasonable  agreement  between  such  theoretical  predictions  and  experimental  data. 
With  respect  to  these  figures,  the  pressurization  of  the  test  chamber  was 
accomplished  by  an  influx  of  nitrogen  gas  which  began  coming  from  the  driver 
chamber  at  1.57  msec.  In  the  actual  vibration  tests  no  propellant  was  burning. 
However,  if  an  M6  propellant  strand  of  cross  sectional  area  0.1186  in.2  had  been 
ignited  at  t  =  1.57  msec  and  had  burned  planarly  in  accordance  with  a  conventional 
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burning  rate  law  (r  =  Bp^ ),  the  theoretically  predicted  transient  test  chamber 
pressures  would  have  been  only  very  slightly  higher  than  those  predicted  pressures 
indicated  in  Figures  1  and  2. 

In  order  to  gain  some  semi-quantitative  and  qualitative  information  about 
the  one-dimensional  mechanical  response  velocity  of  the  exposed  surface  of  a 
propellant  strand  when  subjected  to  forces  caused  by  transient  gas  pressure,  a  finite 
element  analysis  was  applied  to  the  structural  idealization  presented  in  Figure  3.  A 
detailed  description  of  this  structural  idealization  can  be  found  in  Appendix  D.  Two 
aspects  of  the  schematic  in  Figure  3  should,  however,  be  reiterated  here.  First, 
interconnection  nodes  52  through  101  in  conjunction  with  the  attached  springs  were 
intended  to  simulate  the  rigidizing  effect  upon  the  dielectric  material  caused  when  it 
is  bonded  over  its  entire  perimeter  area  to  the  inside  waveguide  walls.  Secondly,  the 
transient  force  F(t)  applied  to  interconnection  node  number  one  was  calculated  by 
taking  the  product  of  cross  sectional  waveguide  area  and  a  time  dependent  test 
chamber  pressure  predicted  by  the  gas  dynamic  model. 

Shown  in  Figures  4-7  are  the  velocities  of  interconnection  node  number 
one  when  subjected  to  four  increasingly  severe  theoretical  test  chamber 
pressurizations.  Each  of  these  pressurizations  was  a  prediction  from  the  gas  dynamic 
model  in  which  the  following  parameters  or  expressions  had  been  fixed,  propellant 
ignition  and  flow  onset  times  (both  at  t  =  0),  test  and  driver  chamber  initial 
volumes  (2.11  in.3  and  88.00  in.3,  respectively),  initial  chamber  temperatures  (both 
540° R),  initial  chamber  gaseous  compositions  (both  100%  nitrogen),  the  explicit 
dependence  of  propellant  burning  rate  on  test  chamber  pressure,  and  strand  cross 
sectional  area  (0.1186  in.2).  Those  parameters  which  were  varied  in  order  to  create 
the  different  pressurizations  were  the  flow  nozzle  throat  area  and  the  initial 
pressures  in  the  test  and  driver  chambers.  Presented  in  Table  1  are  the  identifying 
characteristics  of  the  four  pressurizations  utilized  in  calculating  the  predictions  shown 
in  Figures  4-7. 
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Although  intuitively  obvious,  the  results  in  Figures  4-7  emphasize  that  the 
magnitude  of  the  mechanical  response  velocity  of  the  first  interconnection  node  will 
become  increasingly  large  as  the  severity  of  the  transient  loading  increases.  Even  for 
the  extremely  rigid  structural  idealization  shown  in  Figure  3  and  discussed  in 
Appendix  B,  thv  theoretically  predicted  mechanical  response  velocity  becomes 
roughly  the  same  as  the  propellant  burning  rate  when  the  magnitude  of  the 
pressurization  rate  approaches  that  relevant  to  gun  interior  ballistics.  Further  insight 
into  propellant  strand  mechanical  response  can  be  gained  by  examining  a  highly 
compliant,  though  not  very  realistic,  structural  idealization  of  the  propellant  in  the 
waveguide.  This  idealization  was  constructed  by  eliminating  those  springs  shown  in 
Figure  3  which  were  within  segments  51  through  87.  Application  of  the  same 
transient  loading  associated  with  Figure  4  to  this  soft  system  yielded  the  mechanical 
response  velocity  shown  in  Figure  8.  The  maximum  velocity  indicated  in  Figure  8  is 
167  times  larger  than  the  maximum  velocity  shown  in  Figure  4.  The  information  in 
Figure  8,  as  well  as  that  in  the  preceding  four  figures,  is  based  upon  an  undamped 
response. 


Use  of  the  finite  element  model  in  conjunction  with  the  gas  dynamic 
model  gave  considerable  insight  into  the  magnitude  of  the  mechanical  response  of  a 
solid  propellant  strand  when  subjected  to  a  transient  loading.  This  investigation 
clearly  demonstrated  that  the  mechanical  response  velocity  is  highly  sensitive  to 
strand  confinement  geometry.  Thus,  it  was  apparent  that  it  would  be  a  very 
difficult  task  to  insure  the  sameness  of  mechanical  response  in  the  twin  strand 
differential  system.  As  based  upon  the  theoretical  investigation,  the  physical  situation 
thought  to  give  the  best  chance  for  realizing  the  desired  mechanical  response 
similarity  was  that  wherein  the  complete  perimeter  area  of  each  of  the  two 
propellant  strands  was  epoxy-bonded  to  the  inside  waveguide  walls. 
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3.0  THEORETICAL  MODEL 


A  schematic  of  the  twin  strand  differential  system  along  with  the  relevant 
voltage  waveforms  is  shown  in  Figure  9.  in  expressing  the  voltages  of  the  reflected 
test  and  reference  signals  at  the  monitoring  points,  all  wave  impedance  mismatches 
within  each  arm  except  those  at  the  burning  and  non-burning  surfaces  were  assumed 
to  be  of  negligible  magnitude.  As  depicted  in  this  simplified  figure,  standard  X-band 
waveguides  lead  to  and  go  from  the  propellant  samples  while  the  strands  themselves 
fill  rectangular  waveguides  of  special  cross  sectional  dimension.  Hence,  two  separate 
phase  constants,  and  B,  should  be  used  in  accounting  for  the  accumulative 

W  p 

phase  angle  differences  between  the  source  signal  and  the  two  reflected  signals.  The 
difference  in  phase  angle  between  the  reflected  signals  is 


=  MbT  -  bR)+MCT  "  cr)+WST  "  +  '  5R  >  (I) 

where  <|/ST  and  are  the  phase  shifts  caused  by  a  jump  in  wave  impedance  at 

the  burning  and  non-burning  surfaces,  respectively. 

If  the  lengths  bT ,  bR .  Cj-,  and  cR  are  constant  and  if  4*ST  and  \pSR  do 
not  change  with  time,  then  taking  (d/dt)  of  equation  (1)  yields 

K  =  2VCt  -*r>  (2) 

Assume  that  motion  of  the  burning  surface  results  from  both  a  propellant 
deflagration  and  a  mechanical  or  vibratory  response.  Let  the  non-burning  surface 
undergo  motion  caused  by  only  a  mechanical  response.  Thus, 
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£t 


=  ~  ra+VmT 


(3) 


& 


R 


Vm  R 


(4) 


where  the  burning  rate  ra  is  considered  positiv  -  if  fiT  decreases  during  combustion. 
If  vmR  =  vmT  or  if  both  mechanical  responses  are  zero,  then 


=  -  2/3  r 

a 


(5) 


The  phase  constant  is  given  by 


0 


p 


(6) 


where  Xpp  is  the  wavelength  of  the  microwave  signal  in  the  propellant  filled 
waveguide.  Rearranging  equations  (5)  and  (6)  allows  the  burning  rate  to  be 
expressed  as 


r 


a 


720°  dt 


(7) 


Assume  that  vm  R  and  vmT  are  both  zero  and  let  Hj.  =  Co  at  the  start  of 
burning  while  RT  =  0  at  sample  burnout.  Then  integrating  equation  (2)  over  the 
course  of  a  complete  test  strand  burn  allows  the  total  change  in  relative  phase  angle 
to  be  expressed  as 

I 

►  > 

1 
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A*a  =  W.j-t.j)  =  20p8c 


(8) 


The  phase  wavelength  can  now  be  written  as 


X  =  4?rgo 
Pg  = 


(9) 


Based  upon  microwave  theory ,(28)  the  phase  wavelength 


can  be  expressed  as 


X 


pg 


(10) 


For  the  dominant  TE[0  wave  in  a  dielectric  filled  rectangular  waveguide,  the  cutoff 
frequency  is 


c 


(11) 


If  the  non-burning  surface  is  at  rest,  then  the  previous  development 
reduces  to  a  statement  of  the  Doppler  effect  which  is  the  basis  of  CW  radar.* 2  9) 
The  Doppler  angular  frequency  is  given  by 


coD  =  27T(fR  -  fT)  = 


dt 


4n  dCT 

IT 

pg 


(12) 


» 

t, 
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where  (fR  -  fT)  is  known  as  the  beat  frequency.  In  using  the  conventional  CW  radar 
approach,  the  velocity  of  the  moving  reflecting  surface  is  determined  by  the 
measurement  of  the  beat  frequency.  However,  the  time  interval  required  to 

determine  the  beat  frequency  is  much  too  long  to  have  any  practical  utility  in 
measuring  solid  propellant  burning  velocity  in  a  dynamic  pressure  environment.  For 
example,  if  the  wavelength  of  the  microwave  signal  in  the  propellant  filled 

waveguide  was  1  in.  and  if  the  propellant  burning  rate  was  1  in. /sec,  then  these 
typical  values  would  produce  a  beat  frequency  of  only  2  Hz. 

Shown  in  Figure  10  is  a  phasor  diagram  relevant  to  the  schematic 

presented  in  Figure  9.  However,  the  situation  depicted  in  Figure  10  is  more  general 

since  the  previously  neglected  wave  impedance  mismatches  within  each  branch  of  the 

system  are  now  assumed  to  cause  fixed  reflections  of  significant  magnitude.  Because 
— ^ 

of  the  presence  of  Ax  and  AR ,  the  test  and  reference  signals  are  no  longer 

composed  of  only  the  microwave  reflections  from  the  moving  surfaces,  BT  and  BR . 

— ►  — ► 

Rather,  the  test  and  reference  signals,  CT  and  CR ,  respectively,  must  be  considered 

as  resultant  phasors.  The  phase  angle  difference  between  the  test  -md  reference 

signals  is  \pm  =  (i pm  T  -  R ).  If  it  is  assumed  that  the  fixed  reflections  are  very 

— ►  — ► 

small  compared  to  the  magnitudes  of  BT  and  BR ,  then  i pm  can  be  equated  to 
=  ^aT"  r  )•  It  was  precisely  this  assumption  which  was  used  in  arriving  ai 

equation  (7).  Since  the  experimental  technique  was  intended  to  exploit  equation  (7) 
but  could  only  measure  ^ ,  it  was  very  important  to  try  to  minimize  the 

magnitude  of  the  undesirable  fixed  reflections. 

Shown  in  Figure  1 1  is  a  less  complicated  phasor  diagram  representing  the 
microwave  process  when  a  single  solid  propellant  strand  burns  planarly  within  a 
waveguide.  For  simplicity,  this  phasor  has  been  constructed  by  considering  that  the 
reference  signal  was  of  constant  magnitude  and  that  it  was  rotating  at  the  source 
radian  frequency.  Such  a  reference  signal  could  easily  be  achieved  if  the  reference 
strand  in  Figure  9  was  replaced  with  a  fixed  waveguide  short  that  was  subject  to 
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no  mechanical  excitation.  In  Figure  1 1  the  test  and  reference  signals  are  represented 
— *  — ^ 

by  C  and  V  ,  respectively.  Assuming  that  the  mechanical  response  of  the  propellant 
strand  is  negligible  when  compared  to  its  combustion  response,  the  governing 
equations  relevant  to  this  phasor  diagram  can  be  expressed  as  follows: 


A  sin  a  +  B  sin  ii/  =  C  sin  \i/ 

r  a  ^  rn 

(13) 

A  cos  a  +  B  cos  i b  =  C  cos  i b 

a  x  m 

(14) 

B  =  B0  exp  Oiv!/a) 

(15) 

In  phasor  notion  these  equations  reduce  to 


Ae“ 


+  B  eM  * J 

O 


‘C'  . 


=  Ce'0,  m 


(16) 


As  before,  the  measured  test  microwave  signal  is  the  resultant  of  the  fixed 
background  field  and  the  instantaneous  microwave  field  returning  from  the  moving 
reflective  surface.  The  magnitude  of  this  latter  reflected  wave  depends  upon  the 
following  three  parameters:  p  the  attenuation  constant  for  the  microwaves  in  the 
propellant  fixed  waveguide,  the  instantaneous  value  of  the  relative  phase  angle, 
and  B0  a  constant  associated  with  the  initial  conditions.  The  constant  magnitude  A 
and  the  constant  phase  angle  a  define  the  fixed  microwave  reflection. 

To  determine  the  actual  propellant  regression  rate,  equation  (7)  ought  to 
be  applied  but  only  \pm  can  be  measured  in  the  actual  test  circumstance.  The 
apparent  or  measured  burning  rate  is 


14 


h* 

720°  dt 


(17) 


However,  as  derived  in  Appendix  E,  there  is  a  relationship  between  and  \pm 
which  is 


[A  mA 

1  -  -  cos(^m  -  a)  +  Sin(^m  -  a) 


This  equation  clearly  shows  that  the  measured  burning  rate,  rm ,  exactly  agrees  with 
the  actual  burning  rate,  ra,  only  when  A  =  0. 
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4.0  NONLINEAR  PARAMETER  ESTIMATION 


Interpreting  the  measured  instantaneous  burning  rate,  rm ,  in  the  presence 
of  significant  fixed  reflections  can  be  quite  difficult.  For  the  physical  situation 
represented  by  the  phasor  diagram  in  Figure  11,  a  technique  for  extracting  the 
actual  burning  rate,  ra,  can  be  derived.  Eliminating  C  from  equations  (13-15)  yields 


B,  exp  (pv!' .,)[sin(ii  -  \p  )]  +  A[sin(^  -  a)l  =0  (19) 

O  J  /  7l  u  til 


If  p,  A,  Bo,  and  a  are  at  constant  values  for  a  given  test,  then  for  any  specific 
i i/m  value  equation  (19)  becomes  a  function  of  only  i^a;  that  is. 


F(*,)  =  0 


(20) 


Hence,  for  each  measured  phase  angle.  ,  when  all  other  parameters  have  values, 
there  is  some  actual  phase  angle,  i^a,  which  satisfies  equation  (20).  As  shown  clearly 
in  Figure  12,  equation  (20)  has  multiple  \pd  roots.  Also  indicated  in  this  figure  is 
the  fact  that  for  some  input  parameters  F(J/a)  exhibits  no  zero  crossings  in  certain 
portions  of  the  domain  of  i^a. 

Newton’s  method*30*  can  be  used  to  approximate  the  root  to  within  a 
small  uncertainty.  The  algorithm  terminates  when  the  absolute  value  of  the  change 
in  the  root  approximation  from  one  iteration  to  the  next  falls  below  a  small 
specified  tolerance  e.  To  find  the  correct  root  of  equation  (20)  for  a  given  \pm  and 
a  specific  four  constant  parameter  set,  the  root  finding  procedure  should  be  started 
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with  a  good  first  approximation.  For  regions  close  to  or  at  the  propellant  burnout 
stage,  a  reasonable  assumption  is  that  $a  2;  ,  and  this  should  be  the  first  guess. 

Each  correct  root  must  allow  the  following  physical  constraints  to  be  upheld: 


C 


0 


(21) 


sin(^a  - 
sin (tfm  -  «  ) 
sin(i pa  -  a  ) 


=  same  sign 


(22) 


Note  that  relation  (21)  is  dominant  in  the  constraint  hierarchy.  When  applying 
equations  (22),  which  are  based  upon  the  Law  of  Sines  for  a  plane  triangle,  it  is 
presupposed  that  A,  B.  and  C  are  positive  magnitudes.  When  any  of  the  angle 
differences  in  equations  (22)  are  small  (say  less  than  e).  then  certain  knowledge  that 
equations  (22)  are  upheld  becomes  lacking.  When  this  particular  situation  occurs,  the 
root  is  fixed  until  new  values  of  allow  increases  in  the  relevant  angle 

differences. 


With  the  i p a  root  established  for  the  given  ,  a  theoretical  value  of  C 
can  be  calculated.  However,  there  is  also  a  measured  value  of  C  for  each  \f/ .  These 
facts  form  the  basis  for  the  optimization  process  to  be  described.  Define  the 
objective  function  to  be  minimized  as 


V 


k 


(23) 


17 


The  nonlinear  parameter  estimation  of  the  constants  n,  A,  B0,  and  a  is  achieved  by 
minimizing  the  degree  of  disparity  between  a  set  of  discrete  measured  C  values  and 
the  corresponding  theoretical  C  values.  When  all  in  the  set  monotonically 
increase  or  decrease,  the  objective  function  depends  only  upon  the  four  parameters 
jr,  A,  B0,  and  a.  Hence,  r\  is  a  multivariable  objective  function  amenable  to 
minimization  by  direct  search  techniques.  One  such  numerical  scheme  is  the  Hooke 
and  Jeeves13 1 ,3 2 *  pattern  search;  it  is  this  approach  which  is  described  briefly 
below. 


A  direct  search  method  (the  term  “direct  search”  having  been  coined  by 
Hooke  and  Jeeves)  requires  only  the  evaluation  of  the  objective  function  at  a 
particular  location.  This  is  in  contrast  to  the  gradient  techniques  which  require  both 
function  and  gradient  values  to  be  found  at  any  position.  In  utilizing  the  Hooke 
and  Jeeves  pattern  search  to  locate  the  optimum  of  rf,  both  the  objective  function 
and  the  system  variables  are,  for  practical  purposes,  considered  unrestricted  even 
though  each  of  the  four  independent  variables  is  constrained  to  be  within  a  given 
search  interval.  Since  the  function  r?  is  very  well  behaved  and  the  inequality 
restrictions  are  set  to  be  very  far  from  the  expected  local  optimum,  the  numerical 
optimization  can  be  considered  essentially  unrestricted  in  the  primary  search  space. 

A  generalized  sequential  procedure  for  implementing  the  Hooke  and  Jeeves 
search  is  given  by  the  following  four  items. 

1.  Select  from  within  the  n-dimensional  search  space  a  set  of  feasible 
values  for  the  independent  variables.  Although  this  is  the  initial  base 
vector,  it  is  also  an  old  base  vector  {x°LD}. 

2.  Sequentially  explore  along  each  of  the  \{  direction  the  locations  which 
are  defined  by 


°L<>  ±5 

b,i  I 


(24) 


►  . 
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Temporarily  move  the  new  base  vector  to  each  position  for  which  the 
objective  function  is  improved.  Once  this  local  search  is  performed  for 
all  n-dimensions  in  turn,  the  exploratory  part  of  any  stage  is  complete 
and  a  new  base  vector  |x^EVV|  will  have  been  established.  If,  after  an 
exploratory  stage,  it  is  found  that  |x^EWJ  =  {x°LD|,  the  fixed  step 
sizes  S(  relevant  to  each  xt  should  be  reduced  and  step  2  should  be 
repeated. 

3.  If  |  x^EW  |  #  |xE)ED  |  ,  then  make  the  pattern  or  extrapolative  move 
given  by 

{\I=2KETKLD|  (25) 

Two  cases  arise  after  such  a  move: 

a.  If  the  objective  function  shows  an  improvement  at  |  xp| ,  then 
redesignate  the  last  |x^EW}  as  |x°ED|  and  go  back  to  step  2 
after  noting  that  now  xp  i  must  be  substituted  for  x°ED  in 
equation  (24). 

b.  If  there  is  no  improvement  in  the  objective  function,  then  return 
to  step  2  after  noting  that  the  next  local  exploration  will  begin 
from  the  last  {x^EW  }  which  is  the  best  point  yet  found.  Hence, 
when  back  at  step  2,  the  current  x^EW  must  be  used  in  place  of 
x°ED  in  equation  (24).  Essentially,  all  past  successes  and  failures 
are  forgotten  and  a  fresh  start  of  the  algorithm  begins. 

4.  The  search  sequence  terminates  after  the  fixed  step  sizes  5(  have  been 
reduced  to  some  acceptably  small  values  and  the  local  explorations 
have  failed  to  find  an  enhancement  of  the  objective  function. 


The  acceleration  in  distance,  which  is  the  prime  advantage  of  the  Hooke  and  Jeeves 
pattern  search,  is  apparent  in  item  three  above.  Thus,  although  the  technique  starts 
cautiously  with  relatively  short  pattern  moves,  the  extrapolated  length  rapidly 
increases  with  repeated  successes. 

After  the  Hooke  and  Jeeves  pattern  search  has  been  used  to  determine  the 
optimized  constants  applicable  to  a  given  situation,  the  denominator  of  the 
right-hand-side  of  equation  (18)  can  be  computed.  In  this  computation,  let  the  C 
value  relevant  to  each  <pm  be  considered  the  theoretical  prediction.  With 
experimentally  measured,  can  be  calculated  from  equation  (18)  and  used  in 
equation  (7)  to  find  r  .  Hence,  for  the  physical  situation  represented  by  Figure  11, 
a  method  has  been  derived  for  extracting  the  actual  burning  rate  from  experimental 
data  biased  by  the  presence  of  unwanted,  microwave  fixed  reflections. 
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5.0  EXPERIMENTAL  APPARATUS  AND  TEST  PROCEDURES 


A  schematic  of  the  overall  experimental  system  is  shown  in  Figure  13.  The 
microwave  source  was  an  amplitude-leveled,  backward-wave  oscillator  (BWO)  set  to 
operate  at  10  GHz  and  stabilized  by  a  phase-locked  synchronizer  to  1  part  in  108 
per  second.  All  equipment  relevant  to  the  experiments  is  identified  in  detail  in 
Appendix  F.  An  RF  power  meter  and  a  digital  frequency  meter  were  used  to  give 
visual  verification  of  the  stabilized  frequency  and  output  power  constancy  of  the 
microwave  source.  A  traveling-wave-tube  microwave  amplifier  (TWT),  delivering 
1  watt  output  for  an  input  of  1  mW  or  less,  was  incorporated  in  the  system  to 
insure  that  sufficient  RF  power  was  always  available  to  conduct  an  experiment.  As 
indicated  in  Figure  13,  the  amplified  microwave  source  signal  was  equally  split  by  a 
3  dB  directional  coupler  and  passed  into  the  test  and  reference  arms.  On  the  source 
side  of  this  coupler,  RG-9A/U  coaxial  cable  was  used  to  connect  the  components. 
Within  each  main  arm  of  the  microwave  circuit,  the  microwave  signal  was 
transmitted  via  standard  X-band  waveguide  except  for  a  4.5  in.  segment  of  semi-rigid 
coaxial  cable  linking  the  E-H  tuner  and  the  combustion  fixture. 

With  the  exception  of  a  variable  phase  shifter  located  in  the  first  part  of 
the  test  arm,  the  two  sides  of  the  microwave  circuit  are  identical.  As  the  microwave 
signal  progresses  through  an  arm  of  the  circuit,  10%  of  its  power  is  diverted  via  a 
10  dB  directional  coupler  into  a  cancellation  loop  containing  a  variable  attenuator 
and  a  variable  phase  shifter.  Before  being  fed  back  into  the  main  arm  via  a  3  dB 
reversed  directional  coupler,  both  the  amplitude  and  phase  of  this  diverted  signal  can 
be  adjusted.  At  the  3-port  circulator,  the  main  signal  is  directed  through  an  E-H 
tuner  and  on  to  the  combustion  fixture  where  it  is  reflected  and  begins  the  return 
trip  to  the  circulator.  After  redirection  by  the  3-port  circulator,  the  reflected  signal 
picks  up  the  signal  from  the  cancellation  loop  and  passes  through  a  variable 
attenuator  which  was  used  to  establish  the  initial  input  power  to  the  harmonic 
frequency  converter.  Ferrite  isolators,  which  permit  microwaves  to  pass  in  only  one 
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direction,  were  used  in  both  main  arms  to  reduce  the  effect  of  undesirable 
reflections  upon  the  RF  signals  of  interest.  The  isolator  located  just  downstream 
from  the  microwave  source  prevented  the  entry  of  potentially  damaging  reflections 
into  the  BWO. 

The  reflected  test  and  reference  RF  signals  are  first  processed  by  two 
Hewlett-Packard  instruments,  the  Model  8411 A  harmonic  frequency  converter  and  the 
Model  8410A  network  analyzer.  These  two  devices  working  in  combination  have  a 
nominal  bandwidth  of  10  kHz  and  convert  the  RF  signals  being  measured  into  two 
278  kHz  sine  waves  which  carry  the  same  relative  amplitude  and  phase  information 
as  the  original  RF  signals.  The  amplitude  relationship  between  the  278  kHz  IF 
signals  was  determined  by  a  Hewlett-Packard  model  8413A  phase-gain  indicator,  a 
plug-in  to  the  network  analyzer  which  also  has  a  10  kHz  nominal  bandwidth.  The 
voltage  output  from  this  relative  gain  indicator  was  proportional  to  the  ratio  of  the 
input  RF  voltages  at  the  test  and  reference  sides  of  the  harmonic  frequency 
converter.  The  dc  output  signal  from  this  instrument  has  a  scale  factor  of 
50  mV/dB  where  the  number  of  decibels  is  given  by  20  log,  0(VT/VR ).  The  relative 
gain  indicator  yields  positive  analog  voltage  and  clockwise  meter  deflection  for  ratios 
of  0  to  +  30  dB  and  negative  voltage  and  counterclockwise  meter  deflection  for 
ratios  of  0  to  -30  dB. 

Tire  reference  channel  output  from  the  network  analyzer  is  a  278  kHz  sine 
wave  with  amplitude  fixed  at  about  2  volts  peak-to-peak.  Tire  278  kHz  IF  signal 
coming  from  the  test  channel  output  of  the  network  analyzer  has  an  amplitude 
which  can  vary  between  0  and  3.5  volts  RMS;  the  specific  amplitude  value  depends 
upon  the  test  channel  RF  input  and  a  preselected  test  channel  gain.  This  test 
channel  gain  is  relative  to  the  reference  channel  and  is  established  by  the  use  of 
precision  digital  control  dials  located  on  the  network  analyzer  front  panel;  they  have 
a  69  dB  range  and  a  l  dB  resolution.  For  the  particular  network  analyzer  used  in 
the  experiments,  a  digital  gain  control  setting  of  25  dB  was  needed  to  make  the 
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meter  on  the  relative  gain  indicator  read  about  0  dB  when  test  and  reference  signals 
of  equal  RF  power  entered  the  harmonic  frequency  converter. 

The  phase  angle  relationship  between  the  two  278  kHz  IF  signals  which 
came  from  the  network  analyzer  was  measured  by  a  Dranetz  Engineering  Labs,  Inc., 
phasemeter.  This  unit  was  specially  modified  by  the  manufacturer  to  have  a 
bandwidth  of  7.32  kHz.  The  analog  output  of  this  phasemeter  is  related  to  phase 
angle  by  the  scale  factor  10  mV/deg  while  the  RMS  phase  noise  on  the  output 
signal  is  about  4  m illidegrees  when  operating  at  full  bandwidth.  At  278  kHz  the 
reference  input  to  the  phasemeter  can  vary  from  1.8  to  2.6  volts  peak-to-peak. 
During  this  condition  of  the  reference  input,  the  278  kHz  test  input  to  the 
phasemeter  should  be  between  17  and  177  millivolts  RMS  for  best  results.  To 
realize  this  amplitude  window  on  the  test  channel  input  to  the  phasemeter,  it  was 
necessary  to  amplify  the  278  kHz  IF  input  signal.  Two  means  of  amplification  were 
available;  either  the  digital  test  channel  gain  control  integral  to  the  network  analyzer 
could  be  used  or  external  amplification  of  the  test  IF  signal  could  be  utilized,  as 
indicated  by  Figure  13,  the  latter  approach  was  taken.  Use  of  an  external  20  dB 
amplifier  allowed  the  tracking  of  a  much  larger  increase  in  test  signal  RF  amplitude 
by  the  relative  gain  indicator  than  was  possible  when  the  internal  gain  control  was 
utilized.  Before  entering  the  test  channel  of  the  phasemeter,  the  output  from  the 
20  dB  amplifier  was  conditioned  by  the  simple  passive  high  pass  filter  shown  in 
Figure  14. 

The  phasemeter  analog  output  ranges  were  ±  1.8  volts  dc  or  0  to  3.6  volts 
dc  depending  upon  manual  selection.  Output  from  the  phasemeter  was  conditioned 
by  a  filter-amplifier  having  16  selectable  lowpass  settings.  This  device  was  fabricated 
in-house  and  its  electrical  schematic  is  shown  in  Figure  15.  The  filter  bandwidth  was 
selected  to  be  compatible  with  the  chosen  sampling  frequency  of  the  analog-to-digital 
(A/D)  data  acquisition  system.  This  system  was  built  in-house  and  was  used  to 
record  the  transient  voltage  outputs  from  the  phasemeter,  the  relative  gain  indicator. 
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and  the  two  piezoelectric  quartz  pressure  transducers  mounted  in  the  combustion 
fixture.  Major  components  of  the  A/D  data  acquisition  system  are  identified  in 
Appendix  F.  On  the  1 2-bit  A/D  channels  the  minimum  time  between  samples  was 
30  psec  while  on  the  15-bit  channel  it  was  10  psec.  This  corresponds  to  maximum 
through-put  rates  of  33.3  kHz  and  100  kHz,  respectively.  The  intervals  between 
sampling  by  the  A/D  system  were  selectable  as  follows:  from  10  psec  to  100  psec 
in  10  psec  steps,  from  100  psec  to  1  msec  in  100  psec  steps,  from  1  msec  to 

10  msec  in  1  msec  steps,  and  from  10  msec  to  90  msec  in  10  msec  steps.  For  the 
experiments,  the  input  voltage  ranges  on  the  12-bit  channels  were  ±2.5  volts  and  0 
to  5  volts  while  on  the  15-bit  channel  the  input  voltage  range  was  0  to  4  volts. 
Consequently,  voltage  resolutions  were  1.22  mV  and  0.122  mV  for  the  12-bit  and 

15-bit  channels,  respectively.  The  maximum  number  of  discrete  samples  which  could 
be  stored  in  memory  was  2048  for  each  of  the  four  A/D  channels. 

After  a  test,  all  discrete  data  was  extracted  from  the  static  random  access 
memory  elements,  passed  through  an  audio  cassette  interface  of  in-house 

constaiction,  and  fed  to  an  input  jack  of  a  cassette-type  tape  recorder  within  which 
was  mounted  on  audio  cassette  tape  of  certified  data  quality.  In  the  cassette 

interface,  frequency  shift  keying  (FSK)  translates  the  high  and  low  voltage  logic 
levels  into  two  different  analog  frequencies  each  of  which  is  maintained  for  a 
certain  number  of  cycles.  If  the  audio  cassette  tape  were  to  be  played  in  a  sound 
system,  a  listener  would  hear  a  continuously  changing  pattern  of  two  different 

tones.  To  put  the  test  data  into  a  more  accessible  format,  the  audio  cassette  tape 
was  replayed  while  the  analog  signal  from  an  output  jack  of  the  tape  recorder  fed 
the  cassette  interface  which  in  turn  provided  logical  binary  voltages  to  a  computer. 
This  computer  was  a  model  PDP-11/45  manufactured  by  Digital  Equipment 

Corporation;  it  was  used  to  construct  a  9-track  digital  tape  upon  which  was  written 
the  discrete  voltages  representing  the  four  channels  of  experimental  data.  This 

9-track  tape  was  utilized  by  a  computer  program  which  ran  on  a  CDC-6700 

computer  and  put  the  data  onto  a  demountable  disk  pack.  Once  on  the  disk,  the 
data  could  be  read  via  conventional  FORTRAN  read  statements  and  used  in  further 
data  reduction  procedures. 
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As  indicated  by  Figure  13,  the  four  channels  of  experimental  data  were 
also  recorded  on  a  direct-writing  analog  recorder  having  a  bandwidth  of  125  Hz. 
This  device  was  used  as  a  precautionary  backup  for  the  A/D  system,  to  record 
events  past  the  maximum  sampling  duration  of  the  A/D  system,  and  to  yield  data 
with  which  calculations  could  be  made  of  the  initial  pressures  in  the  combustion 
Fixture  chambers.  This  data,  which  was  the  analog  output  of  the  two  charge 
amplifiers,  resulted  from  a  sudden  combustion  Fixture  depressurization  and  the 

associated  stress  relaxation  in  the  quartz  crystals  mounted  in  the  pressure  gages.  The 

depressurization  to  atmospheric  conditions  was  initiated  immediately  after  a  steady 
or  dynamic  test.  Given  the  difference  between  the  initial  grounded  or  zero  voltage 
and  the  negative  voltage  at  pressure  equilibrium  with  the  atmosphere,  the  pressure 
on  a  gage  at  test  start  could  be  easily  computed. 

To  initiate  sampling  by  the  A/D  system  and  to  trigger  the  capacitive 
discharges  causing  propellant  ignition  and  activation  of  the  valve  inside  the 

combustion  fixture,  a  time  delay  sequencer  was  utilized.  This  instrument  was 
fabricated  in-house  and  its  electrical  schematic  is  shown  in  Figure  16.  The  time 
delay  sequencer  controlled  the  time  of  occurrence  and  sequence  of  three  events.  One 
of  five  selectable  time  resolutions,  the  minimum  being  1  psec,  could  be 

independently  established  for  each  of  two  digitally  selected  delay  periods.  The  two 
channel  capacitive  discharge  circuit  of  in-house  construction  is  shown  in  schematic 
form  in  Figure  17. 

The  combustion  Fixture  shown  in  Figure  18  is  comprised  of  two 
independently  pressurizeable  chambers  separated  by  any  one  of  three  flow  nozzles 
having  throat  diameters  of  0.05,  0.10,  and  0.15  in.,  respectively.  A  flow  nozzle 
throat  is  sealed  by  a  valve  which  is  opened  by  the  action  of  two  m75-A  pressure 
cartridges  manufactured  by  Teledyne  McCormick  Selph  Company.  The  test  chamber 
contains  the  dual  waveguide  which  holds  the  twin  propellant  strands.  A  given 
transient  pressurization  is  caused  by  the  flow  of  nitrogen  gas  from  the  large  driver 
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chamber.  The  time  of  flow  onset  from  the  driver  chamber  as  well  as  hot  wire 
initiation  for  propellant  ignition  can  be  carefully  controlled  by  use  of  the  time 
delay  sequencer.  As  a  safety  precaution  against  an  overpressure  situation,  a 
rupture-disc  assembly  was  positioned  in  the  sidewall  of  fhe  test  chamber  as  shown 
in  the  offset  section  view  in  Figure  19.  The  large  cylindrical  combustion  fixture 
housing  and  most  internal  components  were  constructed  from  A1S1  4340  steel. 
Because  of  the  anticipated  problem  with  surface  corrosion,  all  of  these  steel  parts 
were  coated  with  an  electrolytic  deposit  of  nickel. 

A  dual  propellant  waveguide  was  made  of  copper  and  the  opening  in  each 
side  was  4.000  in.  long  and  had  a  rectangular  cross  section  of  0.298  X  0.398  in. 
Three  such  waveguides  were  fabricated.  The  cutoff  frequency  in  an  air  filled 
waveguide  of  this  dimension  is  14.8  GHz.  To  determine  the  dielectric  constant  of 
the  propellant  chosen  for  use  in  the  experiments,  the  shorted  transmission  line 
method  was  utilized/ 3 3 ,34)  Under  conditions  of  room  temperature  and  atmospheric 
pressure,  the  dielectric  constant  of  M6  propellant  was  measured  to  be  3.238  at 
10.0  GHz.  When  M6  fills  the  propellant  waveguide  the  cutoff  frequency  is 
8.24  GHz.  By  operating  at  10.0  GHz  only  the  dominant  TH,0  mode  will  propagate 
in  the  waveguide  filled  with  M6  while  no  wave  propagation  will  occur  if  the 

waveguide  is  empty.  The  composition,  excluding  residual  volatiles,  of  the  M6 
propellant  used  in  the  experiments  is  presented  in  Table  2. 

The  dielectric  transition  window  is  a  specialized  coaxial-to-waveguide 
adapter  incorporating  two  unique  pressure  sealing  microwave  feedthroughs.  This 

window,  without  the  feedthroughs  in  their  respective  ports,  is  shown  in  Figures  18 

and  19  in  two  90  degree  opposed  sectional  views.  The  body  of  the  dielectric 

transition  window  was  fabricated  from  17-4  PH  stainless  steel.  A  sectional  view  of  a 
microwave  feedthrough  is  shown  in  Figure  20.  In  the  dielectric  transition  window 
the  blind  rectangular  holes  which  mate  with  the  dual  propellant  waveguide  are  filled 
with  Stycast  2741  epoxide  casting  resin  supplied  by  Fmerson  and  Cuming.  Inc.  For 
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the  most  rigid  formulation  of  this  resin,  the  cured  material  has  an  experimentally 
determined  dielectric  constant  of  3.10  at  10.0  GHz.  The  front  face  of  the  dielectric 
transition  window  along  with  two  extra  microwave  feedthroughs  is  shown  in 
Figure  21.  A  feedthrough  connects  a  resin  filled  waveguide  in  the  dielectric 
transition  window  to  a  4.5  in.  section  of  0.0865  in.-dia.  semi-rigid  coaxial  cable 
which  in  turn  connects  to  standard  X-band  waveguide.  This  semi-rigid  coaxial  cable 
plus  the  attached  connectors  required  to  transition  to  an  X-band  waveguide  flange 
are  shown  in  Figure  22. 

To  determine  the  effectiveness  of  RF  power  transmission  by  the  dielectric 
transition  window,  the  voltage  standing-wave  ratio  (VSWR)  applicable  to  each  side 
was  measured  as  a  function  of  frequency  via  a  slotted  section  technique.1 35  *  To 
utilize  this  method  it  was  first  necessary  to  construct  a  microwave  load  which 
would  absorb  all,  or  essentially  all,  of  the  RF  energy  which  exited  the  front  face  of 
the  dielectric  transition  window.  This  special  twin  load,  which  is  shown  in  section 
view  in  Figure  23,  was  fabricated  from  one  of  the  three  dual  propellant  waveguides. 
The  pyramidal  tapers  which  accomplished  the  energy  absorbtion  were  constructed 
from  Eccosorb  MF-124  magnetically  loaded  epoxy  manufactured  by  Emerson  and 
Cuming,  Inc.  The  view  in  Figure  24  shows  the  load  sitting  upon  the  dielectric 
transition  window.  With  the  items  in  this  configuration  and  by  the  use  of  a  model 
415E  Hewlett-Packard  SWR  meter,  measurements  were  made  of  the  VSWR’s  relevant 
to  the  test  and  reference  sides  of  the  dielectric  transition  window.  The  results  of 
these  measurements  are  shown  in  Figure  25.  At  the  operating  frequency  of 
10.0  GHz,  the  test  and  reference  side  VSWR’s  were  1.65  and  1.83.  respectively. 
These  VSWR’s  translate,  respectively,  into  percentages  of  power  transmitted  equal  to 
about  94%  and  somewhat  more  than  91%.  Clearly,  the  dielectric  transition  window 
was  the  major  source  of  the  microwave  fixed  reflections  caused  by  impedance 
mismatches. 

The  dielectric  transition  window  just  discussed  was  the  second  and  not  the 
first  version  which  was  built.  The  first  transition  section  was  a  waveguide-to-waveguide  type. 


Tliis  device  was  a  single  unit  in  which  two  side-by-side  waveguides  changed,  over  the 
course  of  35  in.,  from  standard  X-band  openings  to  special  waveguide  cross  sections 
identical  to  those  in  the  dual  propellant  waveguide.  A  view  of  this  transition  section 
with  an  attached  X-band  waveguide  flange  assembly  is  shown  in  Figure  26.  Although 
considerable  effort  was  expended  to  make  this  waveguide-to-waveguide  transition 
section  capable  of  containing  the  high  pressure  gas  inside  the  combustion  fixture,  no 
success  was  achieved.  Consequently,  this  first  model  of  the  dielectric  transition 

section  was  never  used  in  any  of  the  experiments.  The  complete  failure  of  the  first 
model  of  the  dielectric  transition  window  forced  the  consideration  of  the 
eoaxial-to-waveguide  design  implemented  successfully  in  the  second  model. 

Pressure  sealing  for  a  microwave  feedthrough  was  accomplished  by  the  use 
of  inverted  cones  on  the  face  and  a  copper  washer  on  the  external  shoulder. 

Inverted  cone-type  seals  were  also  utilized  by  the  two  electrical  feedthroughs  which 
passed  the  capacitive  discharges  causing  propellant  ignition  and  explosive  valve 

actuation.  At  interfaces  within  the  combustion  fixture,  the  pressure  was  contained 
by  the  use  of  resilient  metal  c-type  seals.  For  successful  use,  these  resilient  metal 
seals  required  that  the  squeezing  surfaces  have  an  RMS  roughness  of  16  microin.  or 

less.  The  presence  of  an  almost  imperceptible  nick  across  a  seal  line  was  enough  to 

prevent  pressure  containment  by  the  combustion  fixture.  To  still  be  able  to  achieve 
some  degree  of  pressurization  if  such  a  defect  went  undetected,  elastomeric  o-rings 
were  put  behind  the  metal  seals  on  the  major  end  caps  of  the  driver  chamber. 

However,  for  all  the  actual  experiments,  this  precaution  would  have  been 

unnecessary  since  the  c-type  metallic  seals  functioned  perfectly.  To  verify  the  sealing 
capability  of  the  entire  combustion  fixture,  the  driver  and  test  chambers  were 

independently  pressured  to  25  and  15  kpsi,  respectively.  As  indicated  by  two 
Bourdon-type  dial  pressure  gages,  no  significant  loss  of  pressure  from  either  chamber 
occured  during  a  24  hour  time  period. 
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Shown  in  Figure  27  is  a  schematic  of  the  pressurizing  system  used  in  the 
experiments.  Located  within  the  compression  cylinder  is  a  moveable  aluminum  piston 
sealed  along  its  circumferential  boundary  so  that  compressed  gas  above  the  piston  is 
prevented  from  leaking  into  the  oil  below  the  piston.  To  fully  displace  this  piston 
from  the  bottom  of  the  compression  cylinder  to  the  top  required  42  minutes  of 
operation  by  the  oil  pump.  Each  incremental  pressurization  stage  was  therefore  equal 
to  this  time  interval.  The  cumulative  duration  of  booster  pump  operation  depended 
upon  the  desired  final  gas  pressure  in  a  combustion  fixture  chamber.  An  overview  of 
all  equipment  external  to  the  large  protective  shelter  is  presented  in  Figure  28.  At 
the  left  of  this  picture,  near  the  nitrogen  bottle,  is  the  gas  booster  pump 
surrounded  by  a  steel,  box-like  barricade.  A  closeup  view  of  most  of  the  waveguide 
circuitry  is  shown  in  Figure  29.  The  inside  of  the  protective  shelter  is  shown  in 
Figure  30;  in  the  center  of  this  picture  is  the  assembled  combustion  fixture 
connected  to  both  the  microwave  circuitry  and  the  remotely  controlled  valves. 

In  the  steady-state  burning  tests,  a  single  propellant  strand  burned  within 
the  combustion  fixture  in  the  large  closed  chamber  which  was  Conned  by  the 
removal  of  the  explosively  actuated  valve,  the  flow  nozzle  and  its  retaining  elements, 
and  the  test  chamber  volume-filler/flow-detector.  The  maximum  volume  of  this 
closed  chamber  was  204  cubic  in.;  although  on  some  steady  burning  tests,  aluminum 
filler  disks  were  used  to  reduce  this  volume  and  decrease  the  pumping  time  required 
to  reach  a  high  initial  pressure.  Since  the  propellant  strand  did  burn  in  a  closed 
chamber,  there  was  some  increase  in  pressure  during  a  so-called  steady-state  test. 
However,  the  large  volume  of  the  combustion  chamber  kept  this  undesirable  effect 
relatively  small.  When  conducting  a  steady  state  burning  rate  test,  the  reference  arm 
of  the  microwave  circuitry  was  terminated  by  a  fixed  X-bainl  waveguide  short.  This 
action  insured  that  the  reference  signal  voltage  vector  would  be  of  constant 
magnitude  and  rotating  at  the  microwave  source  angular  frequency. 


ft 
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For  every  experiment  whether  it  was  steady  or  dynamic,  the  cancellation 
of  fixed  reflections  was  an  important  aspect  of  test  procedure.  In  the  steady  tests, 
only  the  microwave  fixed  reflections  in  the  arm  relevant  to  the  burning  strand  were 
cancelled.  For  the  dynamic  tests,  however,  the  process  utilized  to  cancel  the  fixed 
reflections  was  successively  applied  to  both  arms  of  the  microwave  system.  In  this 
process,  use  was  made  of  the  special  twin  load  which  had  dimensions  identical  to  a 
dual  propellant  waveguide.  When  the  twin  load  was  positioned  in  the  combustion 
fixture,  each  of  its  elements  could  absorb  essentially  all  the  microwave  radiation 
which  would  have  entered  the  base  end  of  a  sample  M6  propellant  strand.  Thus,  all 
reflected  power  reaching  the  two  entrance  ports  on  the  harmonic  frequency 
converter  would  have  been  caused  by  the  undesirable  impedance  mismatches  in  a 

given  arm  of  the  microwave  system.  Within  each  side  of  the  system  two  types  of 
waveguide  hardware  were  available  for  use  in  reducing  the  fixed  reflections.  These 

items  were  the  E-H  tuner  and  the  components  in  the  cancellation  loop;  they  could 

have  been  utilized  either  separately  or  both  together.  However,  for  all  experiments 
the  attenuators  in  the  cancellation  loops  were  set  at  maximum  attenuation  and  only 
the  E-H  tuners  were  used  to  cancel  fixed  reflections. 

When  cancelling  the  fixed  reflections  within  a  specific  side  of  the  system, 
that  arm  was  terminated  by  one  element  of  the  twin  load.  The  residual  reflected 
RF  energy  in  this  arm  was  fed  into  the  so-called  test  side  of  the  harmonic 
frequency  converter  while  the  main  variable  attenuator  was  set  at  minimum 
attenuation.  In  the  other  arm.  a  fixed  X-band  waveguide  short  was  used  as  the 
termination.  The  main  attenuator  in  this  arm  was  adjusted  until  -23  dBm  of  power 
entered  the  so-called  reference  side  of  the  harmonic  frequency  converter.  By 

monitoring  the  meter  deflection  in  the  relative  gain  indicator,  resetting  the  digital 
gain  control  on  the  network  analyzer,  and  adjusting  the  two  micrometers  on  the 

E-H  tuner,  the  voltage  level  of  the  fixed  reflections  in  the  specific  arm  of  the 

system  being  considered  could  be  reduced  to  a  value  45  dB  less  than  the  initial 

reference  signal  voltage.  Implicit  in  this  procedure  is  the  assumption  that  the 

microwave  fixed  reflections  will  remain  cancelled  after  the  twin  load  is  replaced  by 
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a  propellant  filled  dual  waveguide  and  after  the  space  about  this  waveguide  is 
pressurized  to  some  high  initial  pressure.  Generally,  as  will  be  shown  later,  this  is  a 
poor  assumption. 

The  ends  of  each  M6  propellant  strand  were  initially  planar  with  the  end 
surfaces  of  the  dual  waveguide.  The  body  of  each  strand  was  machined  such  that  its 
cross  sectional  dimensions  were  equal  to  0.288  X  0.388  in.  with  a  tolerance  of 
±0.001  in.  Shown  in  Figure  31,  next  to  a  twin  waveguide,  are  two  strands  of  M6 
propellant  which  have  been  operated  upon  by  a  two-lip  end  mill.  The  specific 
propellant  sample  confinement  characteristics  within  a  waveguide  will  be  discussed  in 
detail  in  the  next  section.  When  a  strand  was  to  be  burned  within  a  waveguide,  a 
0.007  in.-dia.  nichrome  wire  was  maintained  against  the  uninhibited  exposed  face  of 
the  propellant  sample.  This  wire  was  oriented  parallel  to  the  longest  cross  sectional 
dimension  and  along  a  centerline.  A  pyrotechnic  paste  (designated  X-225  by  the 
Naval  Weapons  Center  at  China  Lake,  California)  was  then  painted  on  both  the  wire 
and  the  rectangular  propellant  surface  until  a  layer  about  0.01  in.  thick  was 
established.  The  finely  powdered  materials  in  this  paste  were  potassium  perchlorate 
(KC104),  titanium,  and  boron.  The  binder  was  a  synthetic  rubber  (polyisobutylene 
or  Vistanex)  which  had  been  dissolved  in  a  hexane  solvent.  Table  3  gives  the  weight 
percentages  of  the  paste  components.  To  enhance  the  heat  feedback  to  the 
propellant  surface,  a  strip  of  transparent  tape  was  placed  over  the  coated  M6 
surface.  When  a  170  pf  capacitor  at  175  volts  discharges  through  the  0.8  ohm 
nichrome  wire,  the  pyrotechnic  paste  is  heated,  ignites,  and  causes  the  sympathetic 
ignition  of  the  M6  propellant  surface.  The  Joule  heating  of  the  wire  is  completed  in 
less  than  1  msec. 
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6.0  RESULTS  AND  ANALYSIS 


6.1  Phase  Angle  Resolution 


Figures  32-36  characterize  the  analog  phase  noise  associated  with  the  phase 
angle  measuring  system  depicted  in  Figure  13.  In  determining  the  results  shown  in 
these  five  figures,  the  microwave  source  was  always  stabilized  at  10.0  GHz.  By  the 
use  of  the  phasementer  scale  factor  of  1  deg/10  mV,  the  phase  noise  values  were 
computed  from  the  root-mean-square  voltages  measured  by  the  true  RMS  voltmeter 
identified  in  list  C  of  Appendix  F.  With  respect  to  Figure  13,  the  observations  of 
analog  noise  were  made  at  the  output  of  the  filter-amplifier  which  conditioned 
phasemeter  output.  The  dominant  bandwidth  limiting  components  in  the 
measurement  system  were  the  network  analyzer,  phasemeter,  and  filter-amplifier. 
Respectively,  the  3-dB  bandwidths  of  these  components  are  10  kHz,  7.32  kHz,  and 
16  selectable  frequencies.  Consequently,  the  nominal  system  bandwdith  was 
approximated  by  the  following  equation 


sys 


(26) 


where  each  frequency  is  relevant  to  the  3-dB  amplitude  point.  As  shown  in 
Figure  32,  the  RMS  phase  noise  is  a  cyclical  function  of  measured  phase  angle.  The 
control  of  phase  angle  was  achieved  by  use  of  the  variable  phase  shifter  in  the  first 
part  of  the  test  arm.  For  the  results  in  Figure  32,  both  the  test  and  reference 
signals  were  reflections  from  fixed  X-band  waveguide  shorts;  the  power  of  each 


signal  was  -23  dBm  at  the  input  to  the  harmonic  frequency  converter.  Bandwidth 
control  of  the  system  was  accomplished  by  the  use  of  the  filter-amplifier.  The 
maximum  and  minimum  nominal  system  bandwidths  identified  in  Figures  32-36  are 
associated  with  the  extremes  of  the  16  lowpass  settings  on  the  filter-amplifier. 

To  determine  whether  or  not  the  results  in  Figure  32  were  dependent 
upon  the  spectral  purity  and  stability  of  the  microwave  source,  two  other  types  of 
microwave  source  were  used  in  place  of  the  BWO.  First  to  be  considered  was  a 
reflex  klystron  oscillator  which  was  in  a  Polarad  signal  generator.  Model  1108A.  The 
second  type  investigated  needed  no  input  from  an  external  synchronizer;  it  was  a 
microwave  frequency  synthesizer  manufactured  by  Systron-Donner.  This  device  was 
composed  of  a  driver  unit,  Model  1600A  option  02,  and  an  output  unit. 
Model  161 1A  option  3.  Neither  type  of  microwave  source  yielded  results  different 
from  those  in  Figure  32.  Both  the  cyclical  noise  behavior  and  the  relatively  large 
phase  noise  values  were  caused  by  the  network  analyzer. 

If  the  power  spectral  density  of  noise  is  flat  over  the  frequency  range  of 
interest,  the  RMS  noise  should  be  proportional  to  the  square  root  of  bandwidth.(36) 
The  reasonably  good  agreement  between  the  symbols  and  curves  in  Figures  33  and 
34  show  that  the  phase  noise  can  be  considered  white.  White  noise  implies  only  a 
flat  spectrum  and  does  not  specify  a  distribution  of  amplitudes.  However,  it  seems 
quite  likely  that  the  phase  noise  should  have  a  Gaussian  amplitude  distribution.  If 
this  is  the  case,  then  a  measured  RMS  noise  value  would  equal  o,  the  standard 
deviation  of  a  Gaussian  amplitude  distribution.  Hence,  the  probability  that  the  noise 
voltage  will  be  less  than  o  volts  is  0.683.  Phase  noise  was  also  a  function  of  test 
signal  power  and  an  indication  of  this  dependence  is  given  by  the  results  in 
Figures  35  and  36.  Indicated  power  levels  are  relevant  to  the  input  ports  of  the 
harmonic  frequency  converter.  The  conversion  of  phase  noise  to  a  spatial  resolution 
associated  with  the  translation  of  a  propellant  surface  within  a  waveguide  is 


accomplished  by  multiplying  a  given  RMS  phase  noise  value  by  (2/3p)_I,  the 
reciprocal  of  twice  the  propellant  phase  constant. 

As  based  upon  the  selected  nominal  system  bandwidth  and  the  growth  of 
test  signal  power,  all  steady-state  burning  rate  tests  had  an  RMS  phase  noise  of 
about  0.1  degree  or  less.  On  dynamic  tests,  the  initial  measured  phase  angle  was 
established  so  that  the  phase  noise  was  at  a  minimum  (see  Figure  32);  also,  the  test 
and  reference  signal  power  levels  were  both  set  at  -23  dBm.  Depending  upon  the 
speed  of  pressurization  on  a  given  dynamic  test  and  upon  the  sampling  frequency 
chosen  for  the  A/D  system,  the  nominal  system  bandwidth  was  established  at  3.012, 
4.515,  or  5.341  kHz. 

As  mentioned  earlier,  the  major  contributor  to  phase  noise  and,  hence,  the 
primary  impediment  to  the  realization  of  better  phase  angle  resolution  was  the 
network  analyzer.  Recently,  a  modification  to  the  basic  microwave  detected  system 
used  in  this  investigation  has  been  reported. ,37)  In  the  new  technique,  the  network 
analyzer  was  omitted;  relative  phase  between  test  and  reference  microwave  signals 
was  preserved  in  a  down-conversion  process  utilizing  two  double  balanced  RF  mixers 
and  two  RF  sources  synchronized  at  frequencies  separated  by  500  kHz.  Phase  noise 
in  the  modified  system  was  reported  to  be  limited  to  that  caused  by  the 
phasemeter. 


6.2  Steady-State  Burning  Rate  Tests 

For  the  eleven  steady  tests  conducted.  Table  4  shows  the  relationships 
among  side-wall  inhibiting  techniques,  initial  pressure,  and  the  occurrence  of 
non-planar  burning  within  the  waveguide.  The  strand  configuration  called  type  A  was 
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consistent  with  that  identified  in  the  Design  Rationale  section  as  having  the  best 
chance  for  realizing  the  desired  response  similarity  in  the  twin  strand  differential 
system.  However,  type  A  as  well  as  strand  configuration  types  B  and  C  proved  to 

be  unsatisfactory  for  insuring  that  planar  burning  could  be  maintained  over  the 

entire  4  in.  waveguide  length  when  tests  were  conducted  at  pressures  higher  than 
about  1000  psig.  Even  less  desirable  was  the  fact  that  each  type  of  inhibiting 

technique  showed  one  instance  where  non-planar  burning  occurred  immediately  upon 
propellant  ignition.  The  occurrence  of  significant  non-planar  burning  was  indicated  by 
the  highly  erratic  values  of  measured  phase  angle  and  relative  amplitude  and  by  the 
obvious  increase  in  the  normally  small  pressurization  rate.  Corroboration  of  this 
statement  is  contained  in  Figures  37  and  38  where  raw  data  from  the  eleventh 

steady-state  test  is  presented.  In  Figure  37  the  measured  phase  angle  begins  its 
erratic  behavior  at  about  0.97  second.  At  the  same  time  in  Figure  38  the  pressure 
in  the  closed  chamber  begins  a  rapid  rise  because  of  a  significant  increase  in  burning 
area. 


Burning  rates  were  determined  for  eight  of  the  eleven  steady-state  tests; 
the  only  tests  not  considered  were  the  three  in  which  non-planar  burning  began 

immediately  upon  propellant  ignition.  Perhaps  the  best  way  to  show  how  all  the 
data  was  analyzed  is  to  consider  the  eighth  test  as  an  example  case.  Shown  in 

Figure  39  is  a  plot  of  measured  phaseangle  based  upon  the  linear  connection  of 

2048  discrete  data  points  which  had  been  acquired  by  the  A/D  systems  at  the  rate 
of  one  sample  every  2  msec.  The  raw  data  in  this  figure  represents  the  phase  angle 
change  as  the  propellant  strand  burned  its  full  4  in.  The  hot  wire  causing  propellant 
ignition  began  its  Joule  heating  at  0.020  sec.  In  Figure  40  the  cyclical  nature  of 
the  phasemeter  output  has  been  eliminated  by  a  plot  of  unfolded  measured  phase 

angle.  The  relative  amplitude  change  in  decibels,  as  based  upon  output  from  the 
relative  gain  indicator,  is  shown  in  Figure  41.  Note  that  this  relative  amplitude 
change  is  simply  given  by  the  instantaneous  relative  amplitude  Y  minus  the  initial 
relative  amplitude  Yo;  i.e.,  AY  =  Y-Yq.  The  pressure  history  in  the  closed 
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chamber  for  the  eighth  steady-state  test  is  shown  in  raw  data  form  in  Figure  42. 
Contrast  the  smoothness  of  this  plot  with  the  extreme  pressure  excursion  depicted 
in  Figure  38.  The  peak  pressure  in  Figure  42  occurred  at  strand  burnout.  As  would 
be  expected  for  a  planar  burning  sample,  there  is  a  well  defined  pressure  peak 
which  corresponds  to  a  flat  burning  surface  intersecting  the  plane  of  the  dielectric 
transition  window. 

Although  the  information  in  Figure  41  was  informative,  there  was  a  need 
to  know  the  amplitude  of  the  test  channel  in  absolute  and  not  just  relative  units. 
The  following  discussion  explains  how  the  test  channel  amplitude  was  determined  in 
units  of  volts.  At  the  input  to  the  harmonic  frequency  converter,  the  power 
difference  in  dBm  between  the  test  and  reference  sides  can  be  expressed  as 


rT,dBr 
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null 


X  +  Y 


(27) 


The  term  Xnu(|  represents  the  test  channel  gain  in  dB  which  if  internally  imposed 
by  the  network  analyzer  would  make  the  output  from  the  relative  gain  indicator 
equal  to  zero  volts.  For  the  particular  network  analyzer  used  in  the  experiments, 
Xnul,  equaled  24.8  dB.  The  term  X  indicates  the  dial  setting  of  the  digital  test 
channel  gain  control  on  the  network  analyzer.  Since  the  left-hand-side  of 
equation  (27)  is  in  decibels,  the  following  relation  can  be  written 

20  1°8i°(v1r)  =  X"ull"X  +  Y  (28) 

Recall  that  VR  =  (PRZ)1/2  where  Z  is  the  input  impedance  in  ohms  of  the 
harmonic  frequency  converter  (Z  =  50  il)  and  PR  is  the  reference  side  input  power 
in  watts.  The  operative  equation  can  now  be  written  as 
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VT  =  (PrZ)‘/210 
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where  PR  can  be  determined  from  the  equation  PR  dB  =10  logJ0(PR/l  mW). 
The  use  of  equation  (29)  allowed  Figure  43  to  be  constructed.  Clearly, 
equation  (29)  could  be  applied  to  the  steady-state  tests  because  PR  had  a  known 
fixed  value  during  an  experiment.  Since  the  Y  values  in  equation  (29)  are  discrete, 
the  resulting  VT  values  form  set  of  discrete  measured  test  signal  magnitudes.  These 
VT  .  are  identical  to  the  Cm  ;  written  in  equation  (23). 

To  calculate  the  measured,  apparent,  or  uncorrected  burning  rate, 
equation  (17)  was  applied.  A  prerequisite  to  the  use  of  this  equation  is  knowledge 
of  the  phase  wavelength  Xpg.  An  explicit  expression  for  Xpg  is  given  by 
equation  (9);  however,  computations  cannot  be  directly  undertaken  with  this  relation 
because  Ai//a  is  unknown.  A  solution  to  this  problem  is  to  realize  that  if  \pm  never 
backtracks  on  the  Riemann  surface,  then  Ai//m  will  nearly  equal  Ai p3  and  the  phase 
wavelength  can  be  approximated  with  low  error  by 
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For  the  eighth  steady-state  test,  the  value  of  Xpg  from  equation  (30)  was  computed 
to  be  1.240  in.  Hence,  the  constant  of  proportionality'  Xpg/720°  .needed  in 
equation  (17)  is  1.722  X  10-3  in. /deg.  To  get  the  values  of  dt//m/dt  needed  in 
equation  (17),  central  differences  were  applied  to  the  raw  data  shown  in  Figure  40. 
Multiplying  each  di//m/dt  by  the  proportionality  constant  produced  the  plot  shown 
in  Figure  44.  The  gaps  in  the  plot  are  relevant  to  sections  of  phase  angle  which 
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had  been  created  by  linearization;  this  artifice  was  necessitated  by  the  cyclical 
phasemeter  transitions  from  360°  to  0°.  Because  these  phase  angle  sections  were  not 
“true”  data  they  were  not  considered  in  the  construction  of  Figure  44  or  in  any 
other  analysis  relating  to  phase  angle  derivatives. 

Within  the  same  computer  program  which  constructed  plots  of  raw  data, 
the  distributions  over  frequency  of  both  amplitude  and  energy  were  computed  for 
each  sample  set  of  discrete  physical  data  acquired  during  a  given  test.  The 

techniques  utilized  in  determining  this  spectral  information  are  given  in  detail  in 
Appendix  G.  For  the  phase  angle  data  in  Figure  40,  the  amplitude  and  energy 
spectra  out  to  the  Nyquist  frequency  are  shown,  respectively,  in  Figures  46  and  47. 
Equation  (G-7)  was  the  basis  for  Figure  46  while  Figure  47  was  based  upon 
equation  (G-8).  Analogous  spectral  information  for  the  test  signal  amplitude  data  of 
Figure  43  is  presented  in  Figures  48  and  49.  The  ordinate  scales  in  Figures  46  and 
48  are  in  units  of  decibels  below  the  maximum  spectral  amplitude.  Figures  46-49 
show  very  little  spectral  content  above  30  Hz.  Hence,  for  both  the  phase  angle  data 
and  the  test  signal  amplitude  data  acquired  on  the  eighth  steady-state  test,  a  pass 
band  from  0  to  30  Hz  should  contain  essentially  all  of  the  significant  signal.  This 
statement  was  found  to  also  apply  to  the  pressure  history  relevant  to  the  eighth 
test.  It  must  be  noted  that  the  establishment  of  this  cutoff  frequency  was. 
fundamentally,  a  subjective  decision.  No  technique  of  data  processing  can  by  itself 
allow  one  to  determine  what  is  the  nature  of  the  phenomenon  which  has  been 

sampled.  It  can  only  show  what  is  in  the  sampled  data,  assuming  no  frequency 
aliasing.  Someone  must  decide  what  in  the  sampled  data  is  pertinent  to  the 
phenomenon  and  what  is  not. 

Having  decided  upon  the  appropriate  cutoff  frequency  applicable  to  data 
from  a  given  steady-state  test,  each  sample  set  was  operated  upon  by  a  lowpass 

nonrecursive  digital  filter.  Both  a  smoothing  and  a  differentiating  digital  filter  were 

employed.  Detailed  descriptions  of  these  two  types  of  filters  are  presented  in 
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Appendix  G.  For  the  eighth  steady-state,  the  unfolded  measured  phase  angle, 
combustion  chamber  pressure,  and  test  signal  amplitude  are  shown  in  filtered  form 
in  Figures  50-52,  respectively.  Compare  these  three  figures  with  the  corresponding 
original  raw  data  presented  in  Figures  40,  42,  and  43.  Instead  of  using  conventional 
central  differences  for  di//m/dt,  the  differentiating  digital  filter  was  utilized  to  give 
the  time  derivatives  needed  for  the  computation  of  the  measured  or  apparent 
burning  rate.  For  the  eighth  test,  the  product  of  these  filter  determined  derivatives 
and  the  previously  calculated  proportionality  constant  is  shown  in  Figure  53.  In 
comparing  this  figure  with  Figure  44,  note  that  the  burning  rates  based  upon  central 
differences  are  indicated  up  to  the  essential  burnout  time  of  3.398  seconds  while 

those  based  upon  the  digital  filter  are  not  presented  after  3.328  seconds. 

Even  though  the  plot  in  Figure  53  is  less  ''noisy”  than  that  in  Figure  44. 

the  apparent  burning  rate  still  has  the  significant  undulations  caused  by  the  presence 

of  fixed  microwave  reflections.  The  existence  of  these  fixed  reflections  are  perhaps 

better  appreciated  when  the  results  in  Figure  52  are  interpreted  within  the  context 

— ► 

of  the  phasor  diagram  in  Figure  11.  As  vector  B  grows  in  magnitude  during  its 

— ►  — ► 

rotation  about  the  tip  of  vector  A,  the  resultant  vector  C  will  be  cyclically 

diminished  and  enhanced  in  magnitude.  This  behavior  is  evident  in  Figure  52  where 
the  magnitude  of  vector  C  is  represented  exactly  by  the  experimental  test  signal 
amplitude. 


In  the  determination  of  desired  actual  burning  rates  from  experimental  data 
which  had  been  affected  by  microwave  fixed  reflections,  the  theory'  discussed  in  the 
Nonlinear  Parameter  Estimation  section  was  applied  to  a  given  steady-state  test.  A 
separate  computer  program  was  written  within  which  this  theory  was  implemented. 
The  key  element  in  the  approach  was  the  reduction  in  the  degree  of  disparity 
between  measured  and  theoretically  predicted  test  signal  amplitudes.  Both  of  these 
amplitudes  were  functions  of  measured  phase  angle.  Hence,  for  the  particular  set  of 
discrete  measured  phase  angles  being  considered,  the  optimization  process  focused 
upon  the  minimization  of  r?,  the  objective  function  given  by  equation  (23).  The 
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results  of  objective  function  minimization  for  the  eighth  test  are  presented  in 
Figure  54.  Filtered  experimental  data  was  always  utilized  in  an  optimization  process 
and  in  Figure  54  this  data  is  relative  to  the  time  interval  from  0.040  to 

3.328  seconds.  With  respect  to  equation  (23),  this  time  period  corresponds  to  a 
summation  index  range  of  i  =  21  to  1665.  Another  way  of  demonstrating  the 
degree  of  objective  function  minimization  for  a  given  test  is  exemplified  by  the 
results  in  Figure  55.  This  figure  is  a  restructuring  of  the  results  in  Figure  54  and  it 
should  be  interpreted  with  the  aid  of  Figure  11.  Each  curve  in  Figure  55  is  a 
nondimensional  phasor  plot  representing  the  path  of  the  tip  of  vector  C;  one  curve 
is  theoretically  based  and  the  other  is  constructed  from  experimental  data. 

As  part  of  the  minimization  of  the  objective  function  t?,  the  constants  n, 
A,  Bo,  and  a  were  necessarily  determined.  Hence,  sufficient  information  was  now 
available  with  which  to  evaluate  equation  (18)  according  to  the  directions  given  at 

the  very  end  of  section  4.0.  Note  that  the  computation  of  \pa  from  equation  (18) 

used  the  values  of  m  given  by  the  differentiating  digital  filter.  The  deduced  \pa 

values  plus  the  value  of  Xp?  given  by  equation  (30)  were  then  utilized  in 

equation  (7)  to  compute  the  actual  burning  rate  ra-  This  theoretically  adjusted 

burning  rate  for  the  eighth  test  is  shown  in  Figure  56.  Compare  this  figure  with 

the  uncorrected  or  measured  burning  rate  presented  in  Figure  53.  In  a  slight 
departure  from  the  focus  on  the  eighth  test,  consider  some  results  from  the  fifth 
steady-state  test  which  are  presented  in  Figures  57  and  58  and  serve  as  another 
example  of  the  significant  differences  between  the  uncorrected  and  theoretically 
adjusted  burning  rates.  The  development  of  the  information  in  Figures  57  and  58 
was  analogous,  respectively,  to  that  for  Figures  53  and  56. 

In  the  computer  program  which  calculated  the  actual  burning  rates,  the 
corresponding  actual  phase  angles  were  also  computed.  Given  the  four  optimized 
constants  and  a  specific  4'm  value,  then  each  value  of  the  actual  phase  angle  \pa 
was  a  root  of  equation  (19).  Calculated  actual  phase  angles  for  the  eighth 


experiment  are  shown  in  Figure  59.  Compare  this  figure  to  Figure  50  and  note  the 
reduction  in  undulations  exhibited  by  the  plot  of  actual  phase  angle.  In  Figure  59 
there  are  a  few  places  where  the  value  of  a  i |/a  root  was  held  fixed  for  a  small 
time  interval.  The  necessity  for  creating  these  slight  breaks  in  curve  continuity  has 
already  been  discussed  in  section  4.0.  Note  in  this  figure  that  the  value  of  the 
actual  phase  angle  at  3.328  seconds  was  just  extended  for  the  remaining  time. 
During  the  time  period  from  0.040  to  3.328  seconds  the  actual  phase  changed  by 
2304.7  degrees  while  the  measured  phase  angle  decreased  2290.9  degrees  during  the 
same  time  span.  Because  these  phase  angle  changes  are  almost  the  same  and  since 
both  relate  to  the  almost  complete  combustion  of  the  original  strand  length,  the 
assumption  that  Ai//m  will  nearly  equal  Ai//a  is  upheld  and,  consequently,  the  use  of 
equation  (30)  is  also  sustained. 

Having  fully  discussed  the  eighth  test  as  an  example  case,  the  discussion  of 
results  will  no  longer  be  particularized  to  this  one  steady-state  experiment.  As  based 
upon  tiie  measured  dielectric  constant  for  M6,  the  propellant  waveguide  dimensions, 
and  the  microwave  operating  frequency,  the  theoretical  phase  wavelength  calculated 
by  use  of  equation  (10)  was  1.158  in.  Hence,  the  theoretically  based  value  of  the 
proportionality  constant  in  equation  (7)  was  1.608  X  10“ 3  in. /deg.  Experimental 
results  for  these  quantities  are  shown  in  Table  7.  The  mean  values  from  this  table, 
plus  and  minus  one  standard  deviation,  are  1.176  ±  0.051  in.  and 
(1.633  ±  0.071)  X  10" 3  in. /deg  for  the  phase  wavelength  and  proportionality  constant, 
respectively.  Expressed  as  a  percentage  of  the  mean  values,  the  standard  deviations 
were  about  4.3%.  Note  that  the  mean  experimental  phase  wavelength  is  only  about 
1.5%  higher  than  the  theoretical  value.  Although  the  data  is  sparse,  it  appears  that 
the  experimental  phase  wavelength  was  not  a  function  of  pressure  in  the  range  from 
1  atmosphere  to  about  10,000  psig.  This  experimental  result  in  conjunction  with  the 
theoretical  prediction  from  equation  (10)  seems  to  indicate  that  the  dielectric 
constant  of  the  relatively  stiff  M6  propellant  was  not  significantly  affected  by  strand 
compression  when  under  a  certain  range  of  hydrostatic  pressures. 
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Although  significant  non-planar  burning  did  not  occur  on  the  sixth  test,  no 
results  for  this  experiment  are  presented  in  Table  7.  As  shown  by  the  plot  in 
Figure  60,  the  measured  phase  angle  did  backtrack  upon  the  Riemann  surface  and, 
thus,  Ai^/a  could  not  be  approximated  by  .  This  lack  of  monotonic  behavior  in 
was  caused  by  the  presence  of  large  microwave  fixed  reflections.  Another 
consequence  of  the  large  fixed  reflections  is  shown  in  Figure  61  where  raw  data 
from  the  sixth  test  is  plotted.  Notice  that  for  some  values  of  the  measured  phase 
angle  there  are  multiple  values  of  test  signal  amplitude.  Contrast  Figure  61  with  the 
more  typical  results  in  Figure  54.  For  regions  where  the  values  of  \pm  are  not 
inonotonic  with  respect  to  time,  the  objective  function  defined  in  equation  (23)  is 
inadequate  for  use  in  calculating  an  applicable  set  of  optimized  constants.  The 
objective  function  definition  presumed  that  for  a  given  ^  there  would  be  one  Cm 
and  one  C(.  Hence,  for  the  sixth  test,  equation  (23)  was  applied  only  over  that 
range  of  summation  indices  which  excluded  ambiguity  in  values  of  Cm  and  C(. 

Consider  the  attenuation  of  a  TE)0  wave  in  a  rectangular  waveguide  filled 
with  an  imperfect  dielectric.  The  attenuation  constant  ad  in  units  of  nepers/in.  is 
given  by(28) 


a 


d 


(31) 


where  the  wave  number  k  can  be  written  as  k  =  27rffiue)1^2  and  where  the  ratio  of 
complex  and  real  permittivities  can  be  interpreted  as  the  loss  tangent;  i.e., 
e"/e'  =  tan  6.  If  it  is  assumed  that  n  =  no  and  if  it  is  recognized  that  e  =  e{  eo 
and  c  =  (j u0e0)-  1/2 ,  then  equation  (31)  can  be  recast  in  the  following  form 


► 
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7rf  (er)*  ^ 2  tan  6 

f wr 


(32) 


For  the  specific  situation  of  the  microwave  experiments  with  M6  propellant, 
er  =  3.238,  f,  =  8.240  GHz,  and  f  =  10.000  GHz.  Thus,  for  this  particular  case 
equation  (32)  can  be  rearranged  to  yield 


tan  6  = 


ad 

8.453 


(33) 


when  ad  is  in  the  units  of  nepers/in. 

Consider  equation  (15)  wherein  the  value  of  \pa  governs  the  magnitude  of 

— ► 

B,  the  vector  representing  the  moving  surface  microwave  reflection.  Suppose  instead 
that  this  magnitude  dependence  is  explicitly  based  upon  the  instantaneous  length  ^ 
of  the  test  strand;  i.e.,  B  =  Boexp(p'2T)  where  y  is  a  new  attenuation  constant.  If 
complete  strand  combustion  is  assumed,  then  p'2o  =  pA\//a  and  y'  =  47rju/Xpf,.  While 
in  the  propellant  waveguide,  the  TEJ0  wave  is  attenuated  both  on  its  way  to  and 
after  its  reflection  from  the  moving  propellant  surface.  Thus,  y!  =  -  2ad  because  ad 
in  equation  (31)  relates  only  to  a  one-way  loss.  An  expression  for  the  original 
attenuation  constant  y  can  now  be  written  as 


M  = 


ad^-pg 

2n 


(34) 


For  the  case  of  complete  strand  combustion,  assume  that  the  total  increase 
in  measured  test  signal  amplitude  approximately  equaled  the  total  increase  in  the 
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magnitude  of  vector  B.  This  is  a  reasonable  assumption  and  it  allows  values  of  <*d 
to  be  computed  directly  from  the  total  decibel  increase  in  test  signal  amplitude;  i.e., 
2ad  —  ( AdB/CQ )( 1  neper/8.686  dB).  With  ad  known,  equations  (33)  and  (34)  can 
be  evaluated.  Using  this  approach  values  of  tan  5  and  p  were  computed  for  the 
steady-state  tests;  the  results  are  shown  in  Table  8.  Since  individual  experimental 
values  of  Apg  could  not  be  determined  for  tests  4,  5,  6.  7,  10,  aid  11,  the 
calculation  of  ju  relevant  to  each  of  these  tests  was  based  upon  the  mean  of  those 
experimentally  determined  phase  wavelengths  given  in  Table  7.  The  mean  of  all 
values  of  p  in  Table  8  is  -  0.0516  nepers/rad  and  the  associated  percent  standard 
deviation  is  14%.  An  indication  of  the  invariant  nature  of  the  microwave  attenuation 
in  an  M6  propellant  strand  is  given  by  observing  the  variation  in  total  test  signal  dB 

change  for  each  steady-state  test  except  number  one.  For  these  ten  tests,  the  mean 

value  of  the  total  dB  change  was  19.27  dB  and  the  standard  deviation  expressed  as 
a  percentage  of  the  mean  was  16%.  The  loss  tangent  corresponding  to  this  mean  dB 
change  was  0.0328.  Note  that  although  the  propellant  samples  were  stored  in  an  air 
conditioned  room,  no  additional  attempt  was  made  to  control  the  environmental 
humidity.  Since  the  M6  propellant  does  exhibit  some  hygroscupicity,  the  propellant 

samples  used  in  the  tests  may  have  had  slightly  differing  amounts  of  retained  water 

vapor  and  this  might  have  influenced  the  microwave  attenuation  in  the  strands. 

Consider  now  the  results  from  the  nonlinear  parameter  estimation  theory  as 
applied  to  the  steady-state  experiments.  The  optimized  constants  found  by  the  use 
of  the  Hooke  and  Jeeves  pattern  search  are  shown  in  Table  9.  Notice  that  the 
values  of  Bo  in  this  table  are,  generally  not  much  greater  than  the  values  of  A. 
Since  the  initial  value  of  i p3  was  positive,  equation  (15)  shows  that  the  magnitude 
of  B  at  the  start  of  strand  burning  will  be  even  less  than  Bo.  Only  for  test  number 
eleven  was  the  initial  B  value  more  than  ten  times  larger  than  the  value  of  A. 
Generally,  the  microwave  fixed  reflections  which  had  been  cancelled  during  the  early 
stage  of  an  experiment  showed  a  significant  magnitude  growth  by  the  time  the 
propellant  strand  was  ignited.  The  inability  to  keep  the  microwave  fixed  reflections 
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cancelled  is  attributed,  primarily,  to  the  behavior  of  the  coaxial  waveguide 
components  which  linked  the  standard  X-band  waveguide  to  the  specially  sized 
rectangular  waveguide  in  the  dielectric  transition  window.  That  is,  these  components 
have  interfaces  which  are  susceptible  to  impedance  changes  induced  by  the  relative 
motion  caused  by  increased  stress,  temperature  variation,  or  slight  vibration. 

Particularly  sensitive  interfaces  are  thought  to  have  been  the  exit  plane  of  a 
microwave  feedthrough  and  the  junction  formed  when  the  connector  on  the 
semi-rigid  coaxial  cable  was  attached  to  the  subminiature  jack  on  a  feedthrough. 

Results  for  test  number  one  are  not  given  in  Table  9  because  the  A/D 

system  did  not  function  properly  on  this  experiment  and  no  discrete  data  was 
acquired.  For  tests  2,  3,  6,  8,  and  9  no  significant  non-planar  burning  occurred,  and 
the  values  of  p  based  upon  optimization  agreed  quite  well  with  those  previously 

calculated  on  the  basis  of  the  total  dB  change  in  test  signal  amplitude.  The  mean 
value  of  the  optimized  attenuation  constants  for  these  five  tests  is  -  0.052447  nepers/rad 
and  the  standard  deviation  expressed  as  a  percentage  of  the  mean  is  1 1%.  Because 
planar  burning  ceased  after  a  short  time  on  tests  5  and  11,  there  was  insufficient 
data  available  to  allow  the  optimization  technique  to  converge  to  p  values  which 
accurately  defined  the  true  nature  of  the  exponential  amplitude  growth  of  the 
moving  surface  microwave  reflection.  Hence,  in  the  optimization  runs  for  each  of 
these  two  tests  the  search  interval  for  p  was  made  extremely  narrow  and  was 

centered  about  the  value  of  p  previously  found  from  the  total  dB  change  in  test 
signal  amplitude.  In  these  two  cases,  the  objective  function  ri  was  essentially  only  a 
function  of  Bo,  A,  and  a. 


In  Table  10  a  comparison  is  made  of  the  mean  steady-state  burning  rates 
and  standard  deviations  for  raw.  filtered,  and  theoretically  adjusted  data.  In 
calculating  the  ir  ‘antaneous  burning  rates  for  tests  5.  6,  and  11,  use  was  made  of 
the  mean  proportionality  constant  given  by  the  results  in  Table  7.  A  summary  of 
the  primary  results  for  the  steady-state  tests  is  presented  in  Table  11.  The  mean 
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burning  rate  for  test  number  one  as  shown  in  Table  1 1  was  actually  unrelated  to 
instantaneous  measurements.  Instead,  this  entry  was  based  upon  total  measured  phase 
angle  change  and  total  time  interval  of  burning.  Figure  62  graphically  displays  the 
most  important  information  in  Table  11.  In  this  figure,  the  lines  through  the 
microwave  results  and  through  the  data  from  Grollman  and  Nelson*45 '  were  each 
based  upon  a  least  squares  linear  regression  analysis.  For  the  analysis  performed  on 
the  eight  microwave  data  points,  the  correlation  coefficient  was  0.9991.  Curve  A  in 
Figure  62  represents  typical  closed  bomb  results  for  M6  propellant  at  90°F.(46) 
Although  not  specifically  stated  by  Grollman  and  Nelson,  it  is  assumed  that  curve  C 
is  relevant  to  results  for  M6  propellant  at  “room  temperature.”  The  microwave 
results,  curve  B,  are  for  M6  propellant  at  70±2°F. 

The  interpretation  of  all  the  information  shown  in  Figure  62  must  be 
done  within  a  framework  of  knowledge  about  the  historical  burning  rate  results 
from  strand  burners  and  closed  bombs.  It  should  be  understood  that  substantial 
differences  exist  between  data  obtained  by  different  investigators  using  either  closed 
bomb  or  strand  burner  techniques.  Quoting  from  a  recent  JANNAF  Burn  Rate 
Measurements  Panel,  “there  is  general  lack  of  agreement  between  strand  burner  and 
closed  bomb  data  as  well  as  between  laboratories  using  closed  bombs  or  strand 
burners.’*47*  Some  but  not  all  of  these  differences  may  be  associated  with 
compositional  variability  in  a  given  propellant.  For  example,  a  difference  in  residual 
solvent  level  ot  one  percent  has  been  noted  to  result  in  approximately  a  ten  percent 
change  in  the  burning  rate.*4 8 1  Fortutitous  or  not,  the  microwave  results  in 
Figure  62  lie  midway  between  curves  A  and  C.  For  this  figure  the  uncertainty  in 
pressure  for  the  microwave  data  is  +  25  psi  or  ±  1  %  of  indicated  pressure,  whichever 
is  larger.  An  indication  of  the  uncertainty  in  microwave  burning  rates  is  given  by 
the  vertical  bars  through  each  data  point  except  that  for  test  number  one.  These 
bars  encompass  the  mean  steady-state  burning  rate  plus  and  minus  one  standard 
deviation.  The  average  of  the  percent  standard  deviations  for  seven  microwave  tests 
was  12 %. 
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Some  concern  prompted  by  theoretical  analyses  has  been  expressed  about 
the  possible  detrimental  effects  which  microwave/plasma  interactions  might  have  upon 
the  accuracy  of  burning  rates  determined  from  the  microwave  Doppler/phase  shift 
technique.1 3 ,4 ’4  9 ,50)  Such  effects  would  involve  transient  microwave  reflections 
caused  by  variations  in  flame  zone  ion  and  electron  concentrations.  These  effects 
would  be  manifested  by  a  time  dependence  in  the  wave  impedance  at  the  burning 
propellant  surface;  it  would  no  longer  be  possible  to  consider  the  term  i//ST  in 
equation  (1)  as  time  invariant.  However,  in  a  few  papers  based  upon  experimental 
investigations,  good  arguments  have  been  made  that  microwave  reflections  from  the 
flame  zone  can  be  discounted.1 2,3 7 ,5  l)  In  references  2  and  37,  the  burnout  of  a 
propellant  sample  was  not  accompanied  by  any  jumps  or  abrupt  shifts  in  phase 
angle  or  test  signal  amplitude.  In  reference  51,  the  amplitude  of  the  reflected 
microwave  signal  did  not  change  when  the  burning  process  transitioned  from  a 
normal  propellant  strand  to  one  of  the  same  formulation  but  with  1/2%  KC1  added. 
Experimental  results  from  the  steady-state  tests  with  M6  propellant  were  analogous 
to  those  just  mentioned.  Consider  Figures  50  and  52  wherein  typical  data  shows 
that  neither  the  measured  phase  angle  or  the  test  signal  amplitude  exhibited  sudden 
changes  during  strand  ignition  or  burnout.  If  reflections  from  the  flame  zone  had 
been  important,  then  at  burnout  an  abrupt  change  in  test  signal  voltage  should  have 
been  observed  when  the  plasma  rapidly  decayed  to  a  neutral  gas.  Hence,  either 
essentially  all  the  microwave  energy  was  reflected  by  the  solid/gas  interface  or,  if 
the  energy  did  not  pass  through,  the  concentration  of  ions  and  electrons  was 
insufficient  to  reflect  the  microwave  signal. 

On  each  steady-state  burning  rate  test  a  propellant  ignition  delay  was 
observed.  That  is,  there  was  a  time  interval  between  hot  wire  initiation  and  the 
onset  of  clearly  defined  M6  propellant  combustion.  These  delays  are  indicated  in 
Table  12. 
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6.3  Dynamic  Experiments 


6.31  Vibration  Tests 


The  ability  of  the  twin  strand  differential  system  to  cancel  propellant 
mechanical  responses  was  evaluated  in  a  series  of  vibration  or  non-burning  dynamic 
tests.  By  the  use  of  How  nozzles  with  throat  diameters  of  0.05  and  0.15  in.,  the 
pressurizations  in  these  vibration  tests  were  imposed  at  two  distinctly  different  rates, 
one  slow  and  one  fast.  Pressurizations  were  applied  to  twin  propellant  strands,  a 
single  strand,  and  the  dielectric  transition  window  alone.  In  vibration  tests  with  twin 
propellant  samples,  the  strand  configurations  which  were  investigated  were  types  A 
and  C  as  described  in  Table  4.  Only  the  type  C  configuration  was  used  in 

pressurizations  of  single  strands.  Presented  in  Table  13  is  a  summary  of  the  more 
important  physical  circumstances  associated  with  a  given  dynamic  test.  Note  that 
dynamic  tests  refer  to  both  the  vibration  tests  and  the  transient  burning  tests. 

Common  to  all  the  dynamic  experiments  was  the  use  in  the  test  chamber  of  the 

volume-filler/fiow  deflector.  The  initial  free  volume  of  the  driver  chamber  on  each 
dynamic  test  was  88  in.3  For  the  fourth,  fifth,  and  twelfth  vibration  experiments, 
the  initial  free  volume  in  the  test  chamber  was  3.06  in.3  On  all  other  dynamic 
tests,  the  test  chamber  free  volume  was  initially  2.11  in.3  Other  physical  quantities 
which  governed  the  pressurization  rate  on  a  given  dynamic  test  are  shown  in 

Table  14. 


Data  for  the  ten  vibration  tests  identified  in  Table  13  will  be  presented  in 
a  standard  format  consisting  of  a  sequence  of  five  figures.  Hence,  as  a  means  of 
explaining  this  format  consider  the  first  vibration  test  as  an  example  case. 
Figures  63-67  quantitatively  define  this  test.  The  curves  in  these  five  figures  as  well 
as  the  plots  for  all  other  dynamic  tests,  are  based  upon  discrete  raw  data  which 
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was  processed  by  digital  signal  analysis,  primarily  lowpass  filtering.  Appendix  G  in 
conjunction  with  Tables  6  and  15  fully  specify  the  digital  analysis  technique  applied 
to  data  from  dynamic  tests.  Table  15  lists  the  digital  lowpass  cutoff  frequencies  for 
the  physical  quantities  relevant  to  each  dynamic  test.  As  indicated  in  Table  13  the 
strand  configuration  on  the  first  vibration  test  was  type  A;  this  configuration  was 
initially  believed  to  have  the  best  chance  for  realizing  the  desired  response  similarity 
in  the  twin  strand  differential  system.  The  pressure  history  of  the  test  chamber 
during  test  number  1VT  is  shown  in  Figure  63.  In  this  experiment  as  in  all  the 
other  vibration  tests,  the  capacitive  discharge  to  the  explosive  valve  was  started  at 
t  =  1  msec.  For  test  number  1VT  the  time  interval  between  pressure  cartridge 
initiation  and  the  onset  of  flow  from  the  nozzle  was  0.70  msec.  Based  upon  all  ten 
vibration  tests,  the  mean  value  for  this  flow  delay  was  0.60  msec  and  the  percent 
standard  deviation  was  1 1 7c.  Shown  in  Figure  64  is  the  pressurization  rate  for  the 
first  vibration  test. 

For  all  dynamic  experiments,  power  levels  to  the  test  and  reference  sides 
of  the  harmonic  frequency  converter  were  both  initially  established  at  -23  dBm. 
Hence,  if  the  mechanical  responses  of  the  two  propellant  strands  used  in  test  1VT 
had  been  equal,  the  plot  of  relative  amplitude  shown  in  Figure  65  would  have 
remained  very  close  to  0  dB  for  the  entire  pressurization.  Clearly,  the  results  in 
Figure  65  are  contradictory  to  this  ideal  situation.  An  even  better  indication  of  the 
success  or  failure  of  mechanical  response  cancellation  is  given  by  the  results  of 
Figure  66  wherein  three  phase  angles  are  plotted.  Each  of  the  three  curves  in  this 
figure  begin  at  flow  onset.  Curve  A  is  the  measured  relative  phase  angle.  Curve  B 
represents  a  theoretical  phase  angle  based  upon  the  hypothetical  combustion  response 
of  an  M6  propellant  strand;  the  burning  process  was  assumed  to  be  planar  and  in 
accordance  with  instantaneous  test  chamber  pressure  and  the  steady-state  microwave 
burning  rate  results.  Given  these  assumptions,  this  theoretical  phase  angle  could  be 
interpreted  as  the  actual  phase  angle  associated  with  the  test  signal;  i.e.,  \pa  T  as 
depicted  in  Figure  10. 


With  the  assumed  M6  burning  beginning  at  t  =  0,  each  discrete  value  of 
the  theoretical  phase  angle  was  calculated  from  the  following  equation 


720° 


=  *m.i+ 


Pg 


(35) 


where. 


(A£r)n  (2t  “  ®0)n 


S  B(p2)?  T 


(36) 


In  equation  (35)  the  term  \pm  {  represents  the  measured  phase  angle  at  t  =  0  and 
Xpg  was  the  mean  phase  wavelength  from  the  results  in  Table  7.  The  coefficient  B 
and  the  exponent  a  had  values  of  1.694  X  10~3  and  0.7441,  respectively.  Curve  C 
in  Figure  66  is  a  phase  angle  constructed  from  the  measured  phase  angle  plus  the 
change  in  theoretical  phase  angle;  that  is. 
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Successful  nullification  of  the  mechanical  responses  in  the  twin  strands  would  yield 
curve  A  as  a  horizontal  straight  line  and  curves  B  and  C  as  coincident  traces. 
Negative  slopes  on  the  curves  are  associated  with  apparent  recession  within  the 
waveguide  while  positive  slopes  indicate  apparent  motion  in  the  opposite  direction. 


Shown  in  Figure  67  is  the  final  set  of  results  which  help  to  define  the 
capability  of  the  twin  strand  differential  system  to  nullify  mechanical  responses. 
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Each  curve  in  this  figure  resulted  from  central  differences  applied  to  a  phase  angle. 
To  compute  velocities  from  these  time  derivatives  of  phase  angle,  the  mean  constant 
of  proportionality  from  Table  7  was  utilized.  In  Figure  67,  trace  B  is  a  velocity 
relevant  to  the  theoretical  phase  angle  calculated  from  steady  burning  rate  results. 
Curve  C  is  a  velocity  based  on  the  phase  angle  constructed  from  the  sum  of 
measured  phase  angle  and  the  change  of  theoretical  phase  angle.  Hence,  curves  B 
and  C  in  Figure  67  follow  directly  from  curves  B  and  C  in  Figure  66.  Had  the 
desired  mechanical  response  cancellation  occurred,  curves  B  and  C  in  Figure  67 
would  have  been  coincident.  Clearly,  this  was  not  the  situation  on  the  first  vibration 
test. 


To  verify  that  the  poor  results  on  test  1VT  were  not  an  aberration,  two 
additional  vibration  tests  with  the  type  A  strand  configuration  were  conducted.  The 
pressurization  rates  on  these  twin  strand  differential  experiments  were  analogous  to 
that  in  test  1VT.  Because  ground  loop  effects  in  the  microwave  signal  processing 
equipment  caused  significant  loss  of  data,  the  results  from  these  two  tests  are  not 
presented  in  graphic  form  in  this  report.  However,  it  was  possible  to  deduce  the 

overall  behavior  on  these  two  additional  tests.  Results  for  the  extra  tests  were  very 
similar  to  those  from  test  1VT;  that  is,  they  were  just  as  poor.  It  should  be  noted 
that  no  epoxy  bond-line  failures  were  observed  on  either  the  additional  two  tests  or 
on  the  first  vibration  experiment.  Different  degrees  of  adhesion  or  bond  strength 

between  the  propellant  samples  and  the  waveguide  walls  probably  accounted  for 

most  of  the  mechanical  response  dissimilarity  when  the  type  A  strand  configuration 

was  employed.  Hence,  in  the  remaining  dynamic  tests  attention  was  focused  upon 
the  mechanical  response  behavior  of  propellant  samples  having  the  type  C  strand 
configuration  as  described  in  Table  4. 

To  determine  the  importance  of  the  mechanical  response  associated  with 
the  dielectric  transition  windown  alone,  the  tests  identified  as  4VT,  5VT.  and  12VT 
were  conducted.  On  experiments  4VT  and  5VT  the  pressurizations  were  slow  while 
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on  test  12VT  a  fast  pressurization  was  imposed.  The  results  from  these  tests  are 
presented  in  the  following  three  sets  of  figures:  4VT  in  Figures  68-72,  5VT  in 
Figures  73-77,  and  12VT  in  Figures  78-82.  Each  set  of  figures  is  cast  in  exactly  the 
same  pattern  as  that  just  described  for  the  first  vibration  test  which  served  as  the 
example  case.  Presented  in  each  set  are  the  time  dependent  values  of  pressure, 
pressurization  rate,  relative  amplitude,  phase  angles,  and  velocities.  Note  that  since 
no  propellant  was  present  in  the  dual  waveguide,  both  the  test  and  reference  signals 
were  microwave  reflections  from  the  free  surfaces  of  the  casting  resin  in  the 
dielectric  transition  window.  For  the  less  severe  pressurization  rate  applied  in  tests 
4VT  and  5VT,  the  mechanical  response  of  the  dielectric  transition  window  was 
significant  but  not  nearly  as  large  as  that  displayed  on  test  12VT  when  the  fast 
pressurization  rate  was  imposed.  To  underscore  this  statement  compare  Figures  72 
and  77  with  Figure  82. 

In  tests  I0VT  and  UVT,  respectively,  slow  and  fast  pressurizations  were 
applied  to  single  M6  propellant  samples  each  of  which  had  been  positioned  within  a 
waveguide  according  to  the  type  C  strand  configuration.  The  reference  signal  in 
these  tests  was  the  microwave  reflection  from  an  X-band  waveguide  short  which  was 
fixed  in  its  location  and  underwent  no  relative  motion  during  a  pressurization. 
Figures  83-87  show  the  results  from  test  10VT  and  Figures  88-92  define  test  1 1VT. 
Particularly  interesting  results  from  these  two  tests  are  contained  in  the  phase  angle 
plots  shown  in  Figures  86  and  91.  Relative  to  the  initial  conditions,  both  measured 
phase  angles  in  these  figures  show  an  increase  in  magnitude  by  the  time  pressure 
equilibrium  begins.  These  phase  angle  increases  indicate  an  apparent  growth  or 
extension  of  the  propellant  strand  within  the  waveguide.  Although  unexpected,  this 
phenomenon  easily  could  have  been  caused  by  some  particular  confinement  geometry 
created  when  the  propellant  strand  responded  to  a  pressurization.  That  is,  the 
physical  situation  about  the  strand  could  have  allowed  the  formation  along  the 
sample  perimeter  area  of  time  dependent,  locally  different  tractions  which  acted  to 
advance  the  microwave  reflecting  surface. 
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Notice  in  Figures  90  and  91,  respectively,  the  sudden  jumps  in  relative 
amplitude  and  measured  phase  angle  which  occurred  upon  the  relaxation  of  the  fast 
pressurization  rate.  As  shown  by  the  comparison  of  Figures  87  and  92,  the  more 
severe  pressurization  caused  a  larger  mechanical  response  in  the  single  M6  propellant 
strand.  As  was  the  case  for  all  pressurizations  involving  the  type  C  strand 
configuration,  the  propellant  samples  used  in  tests  10VT  and  11VT  did  not  lose 
their  positional  integrity  with  respect  to  a  waveguide;  that  is,  no  epoxy  bond-line 
failures  were  observed.  Also  no  separation  of  the  Fonnvar  coating  from  a  propellant 
sample  was  detected  in  any  vibration  test. 

Two  slow  pressurizations  were  applied  to  the  twin  strand  differential 
system  in  which  the  strand  configuration  was  type  C.  Results  for  tests  6VT  and 
7VT  are  shown,  respectively,  in  Figures  93-97  and  Figures  98-102.  Consider  the 
phase  angle  plots  in  Figures  96  and  101;  notice  in  each  figure  that  only  after  the 
beginning  of  pressure  equilibrium  did  the  slope  of  curve  C  approach  that  of 
curve  B.  Note  that  the  measured  phase  angle  plots  in  these  two  figures  were  closer 
to  being  the  desired  horizontal  straight  lines  than  the  measured  phase  angle  plot  in 
Figure  86.  Thus,  even  though  the  results  were  very  poor  for  the  twin  strand 
differential  system  when  exposed  to  a  slow  pressurization,  they  were  still  slightly 

better  than  the  results  for  a  single  strand  responding  to  a  slow  pressurization. 

The  twin  strand  differential  system  with  the  type  C  strand  configuration 

was  also  evaluated  under  the  conditions  of  fast  pressurization.  Two  such  vibration 

tests  were  conducted;  results  from  test  8VT  are  presented  in  Figures  103-107  while 
Figures  108-112  define  test  9VT.  The  most  obvious  results  from  these  two  tests 
were  associated  with  the  relaxation  of  pressurization  rate  during  the  onset  of 
pressure  equilibrium.  This  phenomenon  was  associated  with  sudden  jumps  in  both 
relative  amplitude  and  measured  phase  angle;  Figures  105-106  and  Figures  110-111 
clearly  show  the  severe  excursions  in  tiiese  measured  quantities.  Compare  the 

velocity  plots  in  Figures  107  and  112  with  the  analogous  plots  in  Figure  92  which 
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are  relevant  to  a  single  propellant  sample  experiencing  a  fast  pressurization.  In 
making  this  comparison  note  that  the  twin  strand  differential  system  showed 
mechanical  response  behavior  as  poor  as  that  exhibited  by  the  single  M6  strand. 

In  summary,  the  twin  strand  differential  system  with  either  the  type  A  or 
C  strand  configuration  was  not  effective  in  cancelling  mechanical  responses  in  the 
test  and  reference  propellant  samples.  For  vibration  tests  which  had  the  same 
physical  circumstances  and  about  the  same  pressurization  rate  (i.e.,  4VT-5VT, 
6VT-7VT,  and  8VT-9VT)  the  results  were,  at  best,  qualitatively  and  not 
quantitatively  reproducible.  Finally,  the  dominant  mechanical  responses  were  in  the 
propellant  strands  and  not  the  dielectric  transition  window. 


6.32  Transient  Burning  Tests 


Although  the  results  from  the  vibration  experiments  had  been  poor,  two 
transient  burning  tests  at  different  pressurization  rates  were  conducted  with  the  twin 
strand  differential  system  and  the  type  C  strand  configuration.  Results  for  the  fast 
and  slow  pressurization  experiments  are  shown,  respectively,  in  Figures  113-117  and 
Figures  118-122.  In  observing  these  figures,  note  that  for  each  test  the  A/D  system 
started  at  t  =  0,  the  hot  wire  which  controlled  test  strand  ignition  began  its 
response  at  t  =  1  msec,  the  pressure  cartridges  in  the  valve  were  initiated  at 
t  =  21  msec,  and  gas  from  the  driver  chamber  began  to  flow  through  the  nozzle  at 
about  21.6  msec.  In  the  preparation  of  each  reference  strand  used  in  tests  1TBT 
and  2TBT,  the  Formvar  coating  was  established  over  the  entire  surface  of  the 
propellant  sample.  This  thin  thermal  barrier  was  intended  to  prevent  reference  strand 
ignition  until,  at  least,  the  end  of  the  significant  pressure  transient  caused  by  flow 
from  the  driver  chamber.  After  study  of  the  strip  chart  records  for  relative 
amplitude  and  phase  angle,  it  was  deduced  that  reference  strand  ignition  occurred  at 
about  t  =  180  msec  on  test  1TBT  and  about  t  =  292  msec  on  test  2TBT. 
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The  histories  of  pressure  and  pressurization  rate  for  test  !TBT  are  shown, 
respectively,  in  Figures  113  and  114.  On  this  experiment  the  test  chamber  pressure 
rose  115  psi  during  the  time  interval  t  =  1  to  t  =  21  msec.  Consistent  with  the 
results  from  vibration  tests  with  fast  pressurization,  both  the  relative  amplitude  and 
the  relative  phase  angle  showed  abrupt  changes  during  the  relaxation  in 
pressurization  rate.  In  Figures  115  and  116  at  about  t  =  26  msec  these  jumps  are 
evident.  The  results  in  these  figures  also  suggest  the  presence  of  non-planar  burning, 
i.e.,  burning  along  the  longitudinal  perimeter  area  of  the  propellant  strand.  The  test 
strand  did,  in  fact,  experience  this  type  of  burning.  The  test  propellant  sample  was 
fully  consumed  at  about  t  =  180  msec,  the  same  time  at  which  the  reference  strand 
began  to  burn.  The  total  relative  amplitude  increase  observed  during  test  strand 
burning  on  experiment  1TBT  was  18.^  dB.  All  of  this  relative  amplitude  increase 
occurred  after  the  events  depicted  in  Figure  115. 

The  final  collection  of  results  for  the  fir^t  transient  burning  test  are  the 

velocities  shown  in  Figure  117.  Curve  B  in  this  figure  represents  a  theoretical 

burning  rate  computed  by  the  use  of  the  instantaneous  test  chamber  pressure  in 
conjunction  with  the  microwave  steady-state  burning  rate  correlation  given  in 
Figure  62.  This  hypothetical  combustion  response  was  assumed  to  start  at 
t  =  1  msec,  the  instant  of  hot  wire  initiation  on  the  test  side.  Curve  A  in 

Figure  117  was  calculated  by  the  direct  application  of  equation  (17)  with  Xpp  given 
by  the  mean  of  the  experimental  va’ues  in  Table  7.  Time  derivatives  of  measured 
phase  angle  were  based  upon  central  differences.  Because  of  the  approach  to 

coincidence  exhibited  by  the  traces  in  Figure  117,  it  can  be  argued  that  the 
propellant  had  just  ignited  or  was  close  to  being  ignited  by  t  =  21  msec.  Such  a 
20  msec  ignition  delay  would  not  have  been  substantially  dissimilar  from  the  results 
m  Table  12.  The  certainty  of  unequal  strand  mechanical  responses  in  addition  to 
fe*r  •ade  non-planar  burning  serve  to  completely  invalidate  the  apparent  experimental 
"  i'  lme  rjte  shown  in  Figure  I  17. 
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The  results  for  the  transient  burning  test  conducted  under  the  conditions 
of  slow  pressurization  are  presented  in  Figures  118-122.  The  information  in  these 
five  figures  exactly  parallels  that  just  presented  in  Figures  113-117.  On  test  2TBT 
between  t  =  1  and  t  =  2 1  msec  the  test  chamber  pressure  increased  by  82  psi.  At 
t  =  89  msec,  which  is  near  the  end  of  the  pressure  transient  caused  by  the  influx 
of  gas  from  the  driver  chamber,  the  test  chamber  pressure  from  Figure  118  was 
7142  psi.  The  last  data  point  on  this  figure  is  at  t  =  184.9  msec  and  the  pressure 
was  9070  psi  at  this  time.  As  indicated  by  predictions  from  the  gas  dynamic  model 
described  in  Appendix  C,  this  dramatic  pressure  increase  could  not  possibly  have 
been  caused  by  planar  burning  of  the  test  propellant  strand.  Hence,  burning  along 
the  longitudinal  perimeter  area  of  the  test  sample  was  definitely  present  on 
experiment  2TBT.  In  fact,  by  analysis  of  the  strip  chart  records  for  relative 
amplitude  and  phase  angle,  it  was  deduced  that  the  test  strand  was  fully  consumed 
by  about  t  =  292  msec.  At  this  time  the  test  signal  amplitude  had  increased  by 
19.2  dB,  a  value  fully  consistent  with  those  presented  in  Table  8.  Most  of  this 

increase  in  relative  amplitude  occurred  after  the  events  depicted  in  Figure  1 20. 

Just  as  in  test  ITBT,  the  reference  propellant  strand  in  test  2TBT  ignited, 
purely  by  chance,  at  about  the  same  time  that  the  test  sample  burned  out.  The 
reference  strand  on  test  ITBT  exhibited  intermittent  instances  of  significant 
non-planar  burning  over  the  course  of  combustion.  However,  the  reference  strand  on 
test  2TBT  burned  planarly  during  almost  all  of  its  period  of  combustion.  The  basis 
for  this  statement  is  the  regularity  of  change  exhibited  by  the  measured  phase  angle 
trace  on  the  strip  chart  record.  With  the  exception  of  opposite  slope,  the 
appearance  of  this  trace  was  analogous  to  that  of  the  data  in  Figure  39.  Over  a 
two  second  time  interval  near  the  end  of  reference  strand  burning,  the  strip  chart 

record  showed  that  the  measured  relative  phase  angle  increased  by  1630  degrees. 

Multiplying  (1630  deg/2  sec)  by  1.633  X  10' 3  in. /deg,  the  mean  proportionality 

constant  form  Table  7.  yielded  1.331  in. /sec.  This  operation  was,  essentially,  the 
application  of  equation  (17)  without  the  minus  sign.  During  this  same  2  second 
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time  period,  the  test  chamber  pressure  was  nearly  steady,  decreasing  only  from 
8350  psig  to  8126  psig.  Hence,  the  linear  average  pressure  relevant  to  this  time  span 
was  8238  psig.  The  coordinate  pair  consisting  of  1.331  in./sec  and  8238  psig  form  a 
point  which  lies  very  near  the  linear  regression  line  through  the  microwave  data  in 
Figure  62.  This  result  helps  to  corroborate  the  burning  rate  correlation  found  from 
the  steady-state  microwave  experiments. 

The  velocities  in  Figure  122  were  calculated  via  the  same  techniques  used 
in  the  construction  of  the  information  in  Figure  117.  The  meaning  of  the  apparent 
burning  rate  shown  in  Figure  122  was  totally  destroyed  by  the  presence  of  unequal 
mechanical  responses  in  the  propellant  strands  and  by  the  essentially  immediate 
occurrence  of  non-planar  test  side  burning. 


i 
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7.0  CONCLUSIONS 


The  microwave  Doppler/phase  shift  technique  was  used  successfully  to 
measure  the  steady  burning  rates  of  solid  M6  gun  propellant  at  pressures  from  500 
to  nearly  10,000  psig.  This  M6  data  was  intermediate  to  that  derived  from 
conventional  closed  bomb  and  strand  burner  tests.  The  mean  experimentally 
determined  phase  wavelength  of  the  microwave  signal  in  the  propellant  Tilled 
waveguide  was  within  about  1 .5%  of  the  corresponding  theoretically  predicted  phase 
wavelength.  The  standard  deviation  of  the  experimental  phase  wavelength  was  just 
over  4%  when  expressed  as  a  percentage  of  the  mean  value.  As  based  upon  five 
data  points  over  the  observed  pressure  range,  this  phase  wavelength  showed  no 
discernible  dependence  upon  pressure.  This  information  is  significant  since  it  is  the 
phase  wavelength  which  establishes  the  proportionality  constant  linking  burning  rate 
with  the  time  derivative  of  relative  phase  angle.  Based  upon  the  demonstrated 
pressure  sealing  performance  of  the  combustion  fixture,  particularly  that  of  the 
specially  designed  dielectric  transition  window,  it  should  be  possible  to  use  the 
microwave  technique  for  determining  the  steady  burning  rates  of  other  propellants  at 
higher  pressures  (~  25,000  psi ).  However,  the  proven  inability  to  reliably  maintain 
small  magnitudes  for  the  microwave  fixed  reflections  would  require  that  a  relatively 
complicated  theoretical  adjustment  be  applied  to  the  raw  data  in  order  to  reduce 
the  uncertainty  in  measured  burning  rates. 

lor  the  steady-state  burning  rate  tests,  propellant  recession  was  spatially 
described  via  a  one-dimensional  theoretical  model.  With  the  use  of  four  optimized 
constants  this  model  yielded  predictions  which  compared  favorably  with  experimental 
observations.  This  agreement  indicates  that  an  aggregate  type  of  planar  burning  can 
exist  within  the  propellant  waveguide.  Supporting  this  statement  are  the  relatively 
small  standard  deviations  associated  with  the  steady-state  burning  rate  results.  It  was. 
of  course,  recognized  that  the  physically  burning  Mb  propellant  surface  was  really 
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not  perfectly  planar  in  a  mathematical  sense.  At  any  instant  of  time,  an  actual 
propellant  burning  surface  might  have  irregularities.  That  is,  it  might  possess  a 
characteristic  surface  roughness,  have  local  pits  or  craters,  or  be  at  some  angle  other 
than  90  degrees  from  the  longitudinal  waveguide  axis.  However,  when  using  the 
microwave  technique,  the  only  allowable  interpretation  of  the  state  of  the  microwave 
reflecting  surface  is  that  it  is  perfectly  planar.  Hence,  application  of  the  microwave 
technique  presumes  time  invariance  in  the  spatial  relationship  between  the  actual 
burning  propellant  surface  and  the  effective  surface  which  reflects  the  microwaves. 
With  respect  to  the  combustion  of  M6  propellant,  the  experimental  evidence  suggests 
that  microwave/plasma  interactions  should  not  seriously  compromise  the  basic 
microwave  Doppler/phase  shift  technique.  For  the  steady-state  tests,  the  effects  of 
any  instantaneous  irregularities  in  the  surface  of  the  burning  M6  propellant  were 
collectively  diminished  by  averaging  a  great  many  instantaneous  burning  rates 
together  in  order  to  yield  one  representative  value. 

On  a  transient  burning  test  the  propellant  burns  only  a  very  small  distance 
during  the  time  interval  of  pressurization  and  individual  instantaneous  burning  rates 
become  extremely  important.  Thus,  in  dynamic  burning  tests  propellant  burning 
surface  irregularities  may.  potentially,  have  significant  impact  on  the  values  of  some 
instantaneous  burning  rates,  liven  if  an  “ideal”  twin  strand  differential  system  could 
satisfactorily  isolate  the  combustion  response  during  a  dynamic  pressure  transient,  to 
what  degree  would  burning  surface  irregularities  compromise  the  meaning  of  the 
velocity  of  the  effective  microwave  reflecting  plane?  This  question  becomes  very 
difficult  to  answer  when  it  is  realized  that  the  required  distance  resolution  of  the 
measuring  technique,  say  3.93  X  10  4  in.,  might  be  much  smaller  than  the  local 
imperfections  in  the  burning  propellant  surface.  Also,  the  growth  and  occurrence  of 
these  imperfections  might  be  significantly  time  dependent  over  the  pressurization 
time  interval.  Finally,  different  strands  of  the  same  propellant  might  exhibit  different 
burning  surface  geometry.  Only  after  many  test  replications  with  an  “ideal"  system 
and  after  considerable  statistical  analysis  could  the  question  be  answered.  The 
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preceding  question  relates  to  a  problem  which  is  inherent  when  a  nearly  microscopic 
measurement  technique  is  applied  for  the  purpose  of  quantifying  what  is  really  a 
macroscopic  phenomenon. 

The  microwave  twin  strand  differential  system  built  and  evaluated  in  this 
research  effort  was  unsatisfactory  for  use  in  the  determination  of  burning  rates  of 
solid  propellant  when  under  significant  pressure  transient.  Two  major  problems  in 
the  system  were  immediately  evident.  They  are  as  follows: 

1.  At  even  slow  pressurization  rates  (say,  40  X  103  psi/sec),  the 
mechanical  responses  of  the  twin  propellant  strands  and,  to  a  lesser 
extent,  of  the  dual  components  in  the  dielectric  transition  window  did 
not  exhibit  sufficiently  similar  synchronized  amplitude  characteristics. 

2.  During  propellant  combustion  within  the  test-side  waveguide,  gross 
departures  from  planar  burning  were  unpredictable  in  their  occurrence. 

Governing  the  existence  of  these  two  intimately  related  problems  were  the 
confinement  geometry  of  a  strand  within  a  waveguide  and  the  method  of  burning 
inhibition  used  on  the  longitudinal  perimeter  area  of  a  propellant  sample.  No 
progress  in  using  a  microwave  Doppler/phase  shift  technique  for  dynamic  burning 
rate  determination  will  occur  until  these  two  major  problems  are  overcome. 

Four  other  deficiencies  in  the  presently  configured  microwave  system  have 
been  identified.  These  less  severe  or  second  order  problems  in  dynamic  burning  rate 
determination  are  as  follows: 

I.  Phase  noise  by  the  network  analyzer  limits  the  spatial  resolution  to  a 
value  at  least  20  times  larger  than  that  associated  with  the  phasemeter 
alone. 
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2.  In  a  dynamic  pressure  environment  there  is  no  certainty  that  the 
so-called  fixed  microwave  reflections  actually  are  fixed;  and  if  they 
were  fixed,  their  probable  large  magnitude  would  necessitate  a 
theoretical  adjustment  to  measured  or  apparent  instantaneous  burning 
rate. 

3.  Although  reasonably  good  indirect  evidence  of  its  occurrence  did  exist, 
the  ignition  time  of  the  reference  strand  is  not  directly  indicated. 

4.  With  respect  to  the  time  periods  of  both  slow  and  fast  pressurizations, 
the  ignition  delay  time  for  the  test  strand  is  significantly  variable. 

Work  toward  the  solution  of  these  four  second  order  problems  would  be  warranted 
only  after  the  resolution  of  the  two  previously  identified  first  order  problems. 
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APPENDIX  A 
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THEORETICAL  PREDICTION  BASED  UPON 
A  SIMPLE  GAS  DYNAMIC  MODEL 
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Figure  1.  Comparison  of  the  Prediction  of  the  Gas  Dynamic  Model  with 
Experimental  Data  from  a  Typical  Slow  Pressurization 
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Figure  2.  Comparison  of  the  Prediction  of  the  Gas  Dynamic  Model  with 
Experimental  Data  from  a  Typical  Fast  Pressurization 
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Figure  8.  Mechanical  Response  Velocity  Caused  by  a  Theoretical  Pressurization  and 
Relevant  to  Node  Number  One  in  a  Less  Rigid 
Form  of  the  Structural  Idealization 
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Figure  14.  Electrical  Schematic  of  Simple  Highpass  Filter 


Electrical  Schematic  of  Filter-Amplifier 


Figure  16.  Electrical  Schematic  of  Time  Delay  Sequencer 
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Figure  25.  Voltage  Standing  Wave  Ratio  (VSWR)  for  the  Dielectric 

Transition  Window 
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A,  B,  &  C:  HAND  OPERATED  VALVES  ON  BOOSTER  PUMP 
VI,  V2,  &  V3:  NORMALLY  CLOSED,  SOLENOID-PNEUMATIC. 


Figure  28.  Overview  of  Experimental  Facility 
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Figure  30.  View  Inside  of  the  Protective  Shelter 


Figure  31.  Dual  Propellant  Waveguide  and  Two  M6  Propellant  Strands 


Figure  32.  Phase  Noise  versus  Phase  Angle 
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Figure  35.  Minimum  RMS  Phase  Noise  versus  Test  Signal  Power 
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Figure  40.  Unfolded  Measured  Phase  Angle  on  the  8th  Steady-State  Test 
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Figure  46.  Hybrid  Amplitude  Spectrum  in  dB  for  the  Unfolded  Measured 
Phase  Angle  on  the  8th  Steady-State  Test 
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Figure  49.  Accumulative  Energy  Spectrum  for  the  Test  Signal  Amplitude  on 

the  8th  Steady-State  Test 
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Figure  54.  Test  Signal  Amplitude  versus  Unfolded  Measured  Phase  Angle  for  the  8 
Steady-State  Test:  Optimized  Theory  Compared  to  Filtered 
Experimental  Data 
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Figure  59.  Actual  Phase  Angle  for  the  8th  Steady-State  Test 
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Figure  61.  Test  Signal  Amplitude  versus  Unfolded  Measured  Phase 
Angle  for  the  6th  Steady-State  Test  (Raw  Data) 
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Figure  62.  Burning  Rate  versus  Pressure  for  M6  Propellant 
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Figure  66.  Measured,  Theoretical,  and  Combined  Phase  Angles  for 

Vibration  Test  1VT 


Figure  69.  Test  Chamber  Pressurization  Rate  on  Vibration  Test  4VT 
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Figure  70.  Change  in  Relative  Amplitude  Between  the  Test  and 
Reference  Signals  on  Vibration  Test  4VT 


Figure  71,  Measured,  Theoretical,  and  Combined  Phase  Angles  for 

Vibration  Test  4VT 
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Figure  81.  Measured,  Theoretical,  and  Combined  Phase  Angles  for 
Vibration  Test  12VT 
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Figure  86.  Measured,  Theoretical,  and  Combined  Phase  Angles  for 
Vibration  Test  10VT 
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Figure  91.  Measured,  Theoretical,  and  Combined  Phase  Angles  for 
Vibration  Test  11 VT 
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Figure  92.  Theoretical  and  Pseudo-Experimental  Velocities  for 
Vibration  Test  11VT 
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Figure  97.  Theoretical  and  Pseudo-Experimental  Velocities  for 
Vibration  Test  6VT 
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B:  velocity  based  on  theoretical 
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steady  burning  rate  results 
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Figure  107.  Theoretical  and  Pseudo- Experimental  Velocities  for 
Vibration  Test  8VT 
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Figure  112.  Theoretical  and  Pseudo- Experimental  Velocities  for 
Vibration  Test  9VT 


A:  apparent  experimental  burning  rate 
B:  theoretical  burning  rate  based  on 
steady  burning  rate  results 


Figure  117.  Theoretical  and  Apparent  Experimental  Burning  Rates  for 
Transient  Burning  Test  1TBT 
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Figure  119.  Test  Chamber  Pressurization  Rate  on  Transient  Burning  Test  2TBT 
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Figure  120.  Change  in  Relative  Amplitude  Between  the  Test  and  Reference 
Signals  on  Transient  Burning  Test  2TBT 


Fi*,r«  122.  Theoretical  and  Apparent  Experimental  Burning  Rates  for 
Transient  Burning  Test  2TBT 


Table  1.  Characteristics  of  Theoretical  Pressurizations 
Used  in  Finite  Element  Analysis 
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Relevant 
Figure  Number 

Init.  Dr. 
Ch.  Press, 
(psig) 

Init.  Test 
Ch.  Press, 
(psig) 

Flow  Nzl. 
Throat  Dia. 
(in.) 

Max  dp/dt 
(psi/sec) 

Time  of 
Max  dp/dt 
(msec) 

4,8 

7500 

1500 

0.05 

1.22  X  105 

32.0 

5 

7500 

1500 

0.15 

9.50  X  105 

4.0 

6 

40000 

1 

0.15 

9.00  X  106 

5.3 

7 

40000 

20000 

0.15 

18.30  X  106 

0.8 

Table  2.  Composition  of  M6  Propellant  Used  in  Microwave  Tests 
Constituent  Weight  Percentage 


nitrocellulose  at 

84.04 

13.15%  nitration 

dinitrotoluene 

11.10 

dibutylphthalate 

3.79 

diphenylamine 

1.07 

Table  3.  Composition  of  Pyrotechnic  Paste  Used  in  Microwave  Tests 


Constituent 

Weight  Percentage 

Potassium  perchlorate 

72.3 

Titanium 

14.8 

Boron 

6.9 

Polyisobutylene  or  Vistanex 

6.0 

Hexane 


Sufficient  amount  to  give 
the  desired  mixture  viscosity 
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Table  4.  Relationships  Among  Inhibiting  Techniques, 

Pressure,  and  Non-Planar  Burning 

Steady-State  Initial  Pressure  Occurrence  of  Significant  Inhibiting  Strand 
Test  Number  (psig) _ Non-Planar  Burning _ Coating _ Config. 


1 

479 

no 

Hysol;  EA-9309 

A 

680 

no 

Hysol;  EA-9309 

A 

3 

793 

no 

Hysol:  EA-9309 

A 

4 

3561 

immediately  upon 
propellant  ignition 

Hysol;  EA-9309 

A 

5 

3797 

after  about  0.7  in. 
of  planar  burning 

Hysol;  EA-9309 

A 

6 

3999 

no 

Shell:  Epon  815  & 
V40  curing  agent 

B 

7 

6460 

immediately  upon 
propellant  ignition 

Shell:  Epon  815  & 
V40  curing  agent 

B 

8 

6349 

no 

Fonnvar  solution 

C 

9 

8772 

no 

Formvar  solution 

C 

10 

2295 

immediately  upon 
propellant  ignition 

Formvar  solution 

C 

11 

2283 

after  about  0.5  in. 
of  planar  burning 

Fonnvar  solution 

C 

A.  The  complete  perimeter  area  of  the  strand  was  coated  with  epoxy  as  was  the  entire  inside  wall  area  of 
the  waveguide.  The  strand  was  then  pushed  into  tire  waveguide  and  the  epoxy  was  allowed  to  cure  to  a 
nominal  bond-line  thickness  of  0.005  in. 

B:  The  epoxy  was  painted  on  the  entire  perimeter  area  of  the  propellant  and  pennitted  to  cure.  The  cured 
coating  was  lightly  sanded  until  its  thickness  just  allowed  the  strand  to  fit  within  the  waveguide. 

The  strand  was  then  pushed  into  the  waveguide  to  within  3/8  in.  of  the  waveguide  end  after  all  but  the 
last  3/8  in.  of  strand  had  been  coated  with  a  film  of  Dow  Coming  DC 4  silicone  dielectric  lubricant.  A 
film  of  Hysol  l.poxi-Patch  was  applied  to  the  inside  waveguide  walls  over  the  3/8  in.  length  which 
would  be  closest  to  the  dielectric  transition  window.  The  strand  was  then  pushed  the  remaining 
distance  into  the  waveguide  and  the  small  area  of  epoxy  was  pennitted  to  cure,  thus  insuring  the 
positional  integrity  of  the  strand  during  subsequent  handling. 

C:  The  propellant  was  dipped  in  the  Formvar  solution  and  allowed  to  dry;  then  the  process  was 
sequentially  repeated  until  tire  coating  thickness  just  allowed  the  strand  to  fit  within  the  waveguide. 
(See  second  paragraph  under  item  B  for  remaining  procedures.) 


The  Fonnvar  solution  was  comprised  of  220  ml  of  toluene.  147  ml  of  ethanol,  and  25  gni  of 
polyvinylformal. 
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Table  5.  Digital  Filter  Parameters  Applicable  to  All 
Discrete  Data  Acquired  on  a  Given  Steady-State 
Burning  Test 


Steady-State 


Test  Number 

f,  (Hz) 

f2  (Hz) 

fs  (Hz) 

Ar  (N.D.) 

N  (N.D.) 

2 

2 

4 

50 

0.04 

200 

3 

8 

16 

200 

0.04 

200 

11 

40 

48 

200 

0.04 

200 

5 

35 

55 

500 

0.04 

200 

6 

25 

45 

500 

0.04 

200 

8 

30 

50 

500 

0.04 

200 

9 

40 

60 

500 

0.04 

200 

Dynamic 

Test  Number 

Table  6.  Digital  Filter 
Discrete  Data  Acquiree 

f  (kHz) 

Parameters  Applicable  to  All 
on  a  Given  Dynamic  Test 

XR  (N.D.) 

N  (N.D.) 

1VT 

1 0.0 

0.010 

200 

4VT 

10.0 

0.010 

200 

5VT 

10.0 

0.010 

200 

6VT 

10.0 

0.010 

200 

7VT 

10.0 

0.010 

200 

8VT 

20.0 

0.010 

200 

9VT 

33.0 

0.012 

200 

10VT 

10.0 

0.010 

200 

1 1  VT 

33.3 

0.012 

200 

I2VT 

33.3 

0.012 

200 

1TBT 

33.3 

0.012 

200 

2TBT 

10.0 

0.010 

200 

VT  =  Vibration  Test 
TBT  =  Transient  Burning  Test 
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Table  7.  Experimentally  Determined  Microwave  Parameters 


Phase  Wavelength  for  Proportionality 

Steady-State  M6  —  Filled  Waveguide  Constant,  (Xpg/720°) 

Test  Number  (in.) _  (in./deg) _ 


1 

1.202 

1.669  X  10" 3 

2 

1.103 

1.533  X  10"3 

3 

1.157 

1.606  X  10” 3 

8 

1.240 

1.722  X  10”3 

9 

1.178 

1.636  X  10" 3 

Table  8.  Attenuation  Constants  and  Loss  Tangents  Based 
Upon  Theory  and  the  Total  Change  in  Test  Signal  Amplitude 


Steady-State 

Test  Number 

Test  Signal 
Amp.  lncr. 

<dB) 

«d 

(nepers/in.) 

(in.) 

M 

(nepers/rad) 

tan  6 
(N.D.) 

1 

N.A. 

N.A. 

1.202 

N.A. 

N.A. 

2 

24.32 

0.3500 

1.103 

-0.0614 

0.0414 

3 

18.95 

0.2727 

1.157 

-0.0502 

0.0323 

4 

17.29 

0.2489 

(1.176) 

-0.0466 

0.0294 

5 

17.87 

0.2572 

(1.176) 

-0.0481 

0.0304 

6 

19.80 

0.2849 

(1.176) 

-0.0533 

0.0337 

7 

22.50 

0.3238 

(1.176) 

-0.0606 

0.0383 

8 

13.50 

0.1943 

1.240 

-0.0383 

0.0230 

9 

17.09 

0.2459 

1.178 

-0.0461 

0.0291 

10 

20.00 

0.2878 

(1.176) 

-0.0539 

0.0340 

1 1 

21.40 

0.3080 

(1.176) 

-0.0576 

0.0364 

(1.176)  =  mean 

wavelength  based 

on  tests  1 , 

2,  3,  8,  and 

9. 

1 
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Table  9.  Constants 

from  Optimization  Solutions 

Steady-State 
Test  Number 

H  (nepers/rad) 

B„  (mV) 

A  (mV) 

a  (rad) 

I 

-0.061875 

34.373 

15.941 

0.289941 

3 

-0.053109 

13.399 

1.497 

0.957626 

1 1 

-0.057884 

39.750 

2.000 

3.621250 

5 

-0.048880 

18.517 

12.488 

5.399102 

6 

-0.051250 

53.500 

34.500 

1.866250 

8 

-0.049000 

27.500 

13.250 

4.773185 

9 

-0.047000 

34.375 

14.750 

0.340625 

Table  10.  Comparison  of  Mean  Steady-State  Burning  Rates 
and  Standard  Deviations  for  Raw,  Filtered,  and  Theoretically 

Adjusted  Data 

raw 

data 

filtered 

data 

theor. 

adj.  data 

Steady-State 
Test  Number 

No.  of  Inst. 
Burning  Rates 

r 

(in. /sec) 

0 

(in. /sec) 

r 

(in. /sec) 

0 

(in. /sec) 

r 

(in. /sec) 

0 

(in. /sec) 

1 

N.A. 

0.172 

N.A. 

N.A. 

N.A. 

N.A. 

N.A. 

-> 

687 

0.235 

0.073 

0.235 

0.067 

0.229 

0.017 

3 

1884 

0.257 

0.048 

0.256 

0.041 

0.255 

0.030 

1 1 

180 

0.529 

0.067 

0.530 

0.068 

0.535 

0.057 

5 

360 

0.848 

0.654 

0.848 

0.652 

0.844 

0.121 

6 

1476 

0.790 

0.399 

0.790 

0.397 

0.843 

0.139 

8 

1531 

1.194 

0.299 

1.193 

0.290 

1.196 

0.143 

9 

1 1  18 

1.470 

0.507 

1.470 

0.504 

1.437 

0.168 
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Table  12.  Propellant  Ignition  Delays 


Steady-State 

Test  Number 

Delay  Between  Hot  Wire  Initiation  and 

Onset  of  Clearly  Defined  M6  Combustion 
(msec) 

2 

40 

3 

50 

11 

20 

-  5 

30 

6 

30 

8 

20 

9 

40 

Table  13.  Physical  Circumstances  for  Dynamic  Tests 

Dynamic 
Test  No. 

Pressurization 

Rate,  Qualitative 

Test  Geometry  and  Strand 

Configuration  per  Table  4 

1  VT 

Slow 

Twin  strand  differential;  Type  A 

4VT 

Slow 

Dielectric  transition  window  only 

5VT 

Slow 

Dielectric  transition  window  only 

6VT 

Slow 

Twin  strand  differential;  Type  C 

7VT 

Slow 

Twin  strand  differential;  Type  C 

8VT 

Fast 

Twin  strand  differential;  Type  C 

9VT 

Fast 

Twin  strand  differential;  Type  C 

10VT 

Slow 

Single  strand  vs  fixed  ref.;  Type  C 

1 1  VT 

Fast 

Single  strand  vs  fixed  ref.;  Type  C 

1 2VT 

Fast 

Dielectric  transition  window  only 

ITBT 

Fast 

Twin  strand  differential;  Type  C 

2TBT 

Slow 

Twin  strand  differential;  Type  C 

VT  =  Vibration  Test 
TBT  =  Transient  Burning  Test 
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Table  14.  Pressurization  Rate  Parameters  for  Dynamic  Tests 


Dynamic 

Test  Number 

Flow  Nzl. 

Throat  Diameter 
(in.) 

Initial  Driver 

Ch.  Pressure 
(psig) 

Initial  Test 
Ch.  Pressure 
(psig) 

1  VT 

0.05 

6565 

1514 

4VT 

0.05 

7185 

985 

5VT 

0.05 

7218 

1376 

6VT 

0.05 

7127 

1504 

7VT 

0.05 

7207 

1539 

8VT 

0.15 

7295 

1385 

9VT 

0.15 

7370 

1620 

10VT 

0.05 

7383 

1407 

1 1  VT 

0.15 

7577 

1576 

I2VT 

0.15 

7380 

1452 

ITBT 

0.15 

7385 

1660 

2TBT 

0.05  7290 

VT  =  Vibration  Test 

TBT  =  Transient  Burning  Test 

Table  15.  Digital  Lowpass  Cutoff  Frequencies  for 
Physical  Quantities  Relevant  to  Dynamic  Tests 

1656 

Dynamic 

Test  Number 

Y  m 

Digital  Lowpass  Cutoff  Freq 

P2  Arel 

uencies  (Hz) 

^.h 

^sum 

1  VT 

250 

250 

250 

250 

250 

4VT 

150 

250 

250 

250 

250 

5VT 

250 

250 

250 

250 

250 

6VT 

250 

250 

250 

250 

250 

7VT 

150 

250 

250 

250 

250 

8VT 

1500 

1250 

1500 

1500 

1500 

9VT 

1500 

1250 

1500 

1500 

1500 

I0VT 

250 

300 

300 

300 

300 

1  1  VT 

1000 

1250 

1500 

1500 

1500 

I2VT 

500 

1250 

1500 

1500 

1500 

ITBT 

1250 

1250 

1250 

N.A. 

N.A. 

2TBT 

200 

250 

250 

N.A. 

N.A. 

VT  =  Vibration  Test 
TBT  =  Transient  Burning  Test 
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APPENDIX  C 

GAS  DYNAMIC  MODEL 
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GAS  DYNAMIC  MODEL 


The  simple  gas  dynamic  model  utilized  in  designing  the  combustion  fixture 
is  outlined  below.  Note  that  subscript  “one”  refers  to  the  driver  chamber  and 
subscript  “two”  indicates  the  test  chamber.  The  Heaviside  unit  function,  H(z),  is 
used  to  write  the  governing  equations  so  that  they  apply  to  flow  both  from  and 
into  the  driver  chamber.  Also,  H(z)  is  used  in  the  consideration  of  the  propellant 
ignition  time  and  the  time  of  flow  onset  from  a  chamber.  The  basic  assumptions 
used  in  the  model  development  are  as  follows: 

1.  The  driver  chamber  k  of  constant  volume 

2.  No  heat  loss  occurs  in  either  chamber  or  in  the  flow  nozzle 

3.  Negligible  velocities  exist  in  both  chambers 

4.  The  propellant  burns  planarly  within  the  waveguide  in  the  test 
chamber 

5.  Two  nonreacting  gas  species  are  present,  nitrogen  gas  and  propellant 
combustion  products 

6.  Thermodynamic  gas  mixture  properties  are  based  upon  mass  averages 

7.  The  Nobel-Abel  equation  of  state  is  applicable  in  both  chambers 

The  governing  equations  in  the  gas  dynamic  model  are  as  follows: 


conservation  of  mass: 


1 


d(ni| ) 


=  WH(t  -  t,  L)[H(p2  -  p,)-  H(p,  -  p2)l 


d(m2 ) 


=  WH(t-  t,,L)|H(p,  -  p2)~  H(p2  -  p, )]  +mpH(t-  t|G) 


conservation  of  energy: 


d(~^-  =  WHO  -  t,,L )lh2H(p2  -  p,)-  h,H(p,  -  p2)l 


d(m2L2)-  =  WHO  -  tI  L)[h1  H(p,  -  p2)-h2H(p2-  p,  )1  +  nipEHO  -  t,G)  (C^) 


conservation  of  species  mass; 


d(  in  i  n  -> )  .  . . .  , , 

- 2-  =  WHO  -  t,.  L)(x2  N  H(p2  -  p,)-  x,  N2H(p,  -  p2)l 

dt 


d(m^  =  WHO- t|C)H(t- t,  L)lx2  pH(p2  -  p,)- x,  pH(p,  -  p2)l  (C-6) 


d(m2-'N  2-  =  WHO  -  t,  ,  )|x,  N  H(p,  -  p2)-  x2-N,H(p2  -  p,)l 
dt 


1% 


d(m2  ,P ) 
dt 


riipHU  -  t,G ) 

-  WH(t-  t|G)H(t-  t,L)(x2  pH(p2  -  p,)-  x,  pH(p,  -  p2)] 


(C-8) 


test  chamber  volume  change; 


d(V2) 

dt 


=  rApH(t-  t,G) 


(C-9) 


test  propellant  strand  length  change; 


Mr) 

dt 


rH(t  -  t,G) 


(C-10) 


In  the  above  equations,  the  term  W  represents  the  mass  rate  of  flow  through  the 
How  nozzle  from  either  direction;  it  is  always  considered  a  positive  quantity. 

In  the  auxiliary  equations  given  below,  the  index  i  can  take  on  values  1 
or  2  in  order  to  identify  parameters  relevant  to  the  driver  and  test  chambers, 
respectively. 

mi  =  PiVi  =  mi.p  +  mi.N  2 


1; 


(C-ll) 
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xi,N  2 


mi.N  2 
mi 


Xi.P  = 


=  X 


i.N2cvN 


+  X.  Cv 

I,p  vr 


R,  = 


xi,N  2RN 


+  x  R 

2  >.P  P 


Cy 


N2 


fN2<T> 


c 


fp(T) 


=  C„.  +  R: 


7j  = 


Ri  T s 


(C-12) 


(C-13) 


(C-14) 


(C-15) 


(C-16) 


(C-17) 


(C-18) 


(C-19) 


(C-20) 


(C-2 1 ) 
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nRTu  Kr 
7-I  7  -  1 

'O  'o 


h,  = 


r  =  BPj 


m. 


=  rVr 


(C-22) 


(C-23) 


(C-24) 


(C-25) 


In  equation  (C-22).  nRTo  is  equal  to  the  “force”  or  the  “force  constant” 

of  the  propellant  and  F.  is  equal  effectively  to  the  total  energy  available  from  a 

unit  mass  of  propellant.' 1 3 1  For  M6  propellant  the  respective  values  of  E,  pp,  and 
rjo  are  3.804  X  106  in.-lbf/lbm,  00571  lbm/in.3,  and  29.92  in.3/lbm."4)  For  the 

temperature  range  of  540  to  9000°  R,  the  constant  volume  specific  heat  for  diatomic 
nitrogen  was  computed  from  published  information."5' 

To  determine  properties  relevant  to  propellant  combustion  products,  a 
simple  computer  model  was  constructed  for  determining  the  thermodynamic 

properties  of  an  equilibrium  mixture  of  gases.  The  following  species  were  assumed  to 
be  present  in  the  combustion  products:  CO,  C02,  H20,  OH,  H,  N,.  NO.  02,  O,  N, 
and  H2.  For  each  species,  a  determination  was  made  of  the  number  of  moles,  mole 
fraction,  mass,  and  mass  fraction.  The  technique  utilized  in  computing  these 
quantities  required  the  solution  of  eleven  nonlinear,  simultaneous  alegbraic  equations 
and  was  highly  analogous  to  the  approach  presented  by  Hill  and  Peterson."6'  Four 
equations  resulted  from  the  requirement  of  conservation  of  mass  for  each  of  the 
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elements  carbon,  nitrogen,  oxygen,  and  hydrogen.  The  other  seven  required  relations 
were  based  upon  the  seven  independent  equilibrium  reactions  which  were  selected 
lor  use  in  the  calculation  scheme.  The  temperature  dependent  equilibrium  constants 
associated  with  these  particular  chemical  reactions  were  taken  from  a  tabulation  by 
Hunt.*17’  Specific  heat  data  for  the  eleven  gas  species  were  taken  from  a  NASA 
publication.*18’  The  mixture  properties  were  based  upon  mass  averages  and  were 
calculated  for  a  variety  of  temperature  and  pressure  combinations.  At  each 
thermodynamic  state,  the  calculated  mixture  properties  were  the  molecular  weight, 
gas  constant,  constant  pressure  and  constant  volume  specific  heats,  and  the  ratio  of 
specific  heats. 

For  the  combustion  products  of  M6  propellant  at  500  atmospheres,  the 
model  yielded  data  which  allowed  the  temperature  dependent  constant  volume 
specific  heat  to  be  represented  as 


cv  =  0.26694  +  0.34656  X  10'4T-  0.33733  X  10"8T2  (C-26) 

*  n 


where  T  can  range  from  540  to  5400°R  and  cv  is  in  units  of  (BTU/lbm  -  °R).  At 

p 

the  isochoric  adiabatic  flame  temperature  for  M6  propellant  (T  =  4626°R),  the 
model  predicted  a  value  for  7o  equal  to  1.247.  As  based  upon  predictions  from  the 

equilibrium  model,  the  average  gas  constant  for  the  M6  propellant  gases  was 

815.19  in.-lbf/lbm  -  °R. 

In  calculating  the  mass  rate  of  flow  through  the  nozzle,  the  expressions 
used  for  W  were  based  upon  derivations  which  required  the  use  of  the  perfect  gas 

law.  This  was  done  to  preserve  the  simplicity  of  the  model  even  though  it  does 

introduce  an  inconsistency  since  the  Noble-Able  equation  of  state  is  presumed  to 
apply  in  the  test  and  driver  chambers.  For  the  situation  of  noil-choked  flow,  the 
relation  utilized  for  W  was 
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W  =  C,  Y  — 
d  [1 


(d 


A, 

)4  j  i  /2  I  -  8cPh(Ph  '  Pe^l 


1/2 


(C-27) 


where  the  expansion  factor  is  given  by 


Y  = 


/  1  -  Pi/DO4  \  WJ 
\l-(d,/Dh)4rV_ 


(C-28) 


Equation  (C-27)  is  the  theoretical  adiabatic  relation  for  the  mass  rate  of  How  of  an 
ideal  compressible  fluid  across  section  A,,  but  modified  by  a  discharge 
coefficient."9*  Note  that  the  subscripts  “h”  and  “2”  refer  to  the  high  and  low 
pressure  sides  of  the  nozzle.  In  the  non-choked  situation,  r  =  pc/ph  is  greater  than 
the  critical  pressure  ratio,  rcrj(.  For  choked  or  critical  flow,  r  <  r.rj(  and  the 
expression  used  for  calculating  W  was 


W  =  CdA, 


i  -i 


2*uVhr. 


1/2 


(C-29) 


This  equation,  without  the  modifying  discharge  coefficient,  yields  the  maximum 
theoretical  mass  rate  of  llow  for  the  one-dimensional,  steady,  isentropic  flow  of  a 
perfect  gas.* 20 *  When  applying  equations  (C-27)  and  (C-29),  the  specific  heat  ratio  is 
taken  to  be  the  value  of  7  relevant  to  the  higher  pressure  chamber. 

The  actual  (low  nozzle  was  assumed  to  be  comprised  of  a  frictionless 
converging  nozzle  feeding  a  long  duct  of  constant  cross  section  where  wall  friction 
was  important.  In  order  to  maintain  the  computationally  simple  aspects  of  the  gas 
dynamic  model,  the  practical  effect  of  this  friction  was  accounted  for  by  reducing 
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the  critical  pressure  ratio  from  the  isentropic  flow  value  and  by  slightly  reducing  the 
discharge  coefficient  from  its  typical  value  of  0.98  for  nozzle/2  1 '22)  Hence,  for 
combustion  fixture  design  purposes,  rcrit  was  considered  0.49  and  Cd  was  set  at 
0.93. 


The  fourth-order  Runge-Kuita  technique* 23 *  was  used  to  solve  the  ten 
simultaneous  first-order  ordinary  differential  equations  which  comprise  the  governing 
equations  of  the  gas  dynamic  model.  The  transient  solution  began  from  a  set  of 
initial  conditions  for  die  dependent  variables  and  continued  until  the  end  of  a 
prescribed  time  interval.  The  immediate  output  of  the  transient  solution  was  the 
time  dependent  values  of  the  dependent  variables;  i.e.,  those  quantities  which  appear 
in  parentheses  on  the  left-hand-side  of  equations  (C-l)  through  (C-10).  Further 
reduction  to  other  parameters  of  interest  was  accomplished  by  the  use  of  auxiliary 
equations  (C-l  1  >  through  (C-25).  By  invoking  equation  (5)  in  the  Theoretical  Model 
section,  it  was  possible  to  construct  another  auxiliary  expression  which  related  the 
instantaneous  length  of  the  burning  test  propellant  strand  to  the  instantaneous  actual 
phase  angle  change.  This  very  useful  relationship  is  given  by 


(C-30) 


where  Nt//a  .  is  the  initial  value  of  the  actual  phase  angle. 
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FINITE  ELEMENT  ANALYSIS 
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FINITE  ELEMENT  ANALYSIS 


Presented  in  Figure  3  is  a  schematic  of  the  structural  idealization  used  to 
represent  the  dielectric  material  comprised  of  a  single  propellant  strand  in  a 
waveguide  and  a  contiguous  section  of  the  casting  resin  which  filled  the  blind 
rectangular  holes  in  the  dielectric  transition  window.  The  total  length  of  the 
combination  of  a  virgin  propellant  strand  (length  4.000  in.)  and  the  mating  section 
of  casting  resin  (length  1.066  in.)  was  5.066  in.  This  total  length  was  divided  into 
50  segments.  Hence,  between  interconnection  nodes  1  through  51  there  were  50 
springs  each  having  a  spring  constant  given  by  LA/li  =  k,.  It  was  assumed  that  the 
casting  resin  had  an  elastic  modulus  and  a  Poisson’s  ratio  which  were  equal  to  that 
of  the  propellant.  In  constructing  the  schematic  of  Figure  3.  it  was  assumed  that 
epoxy  bonded  the  entire  perimeter  area  of  both  the  propellant  strand  and  the 
casting  resin  section  to  the  surrounding  waveguide  sidewalls.  Thus,  interconnection 
nodes  52  through  101  were  fixed  in  space  and  connected  by  springs  to  moveable 
nodes  1  through  50.  This  aspect  of  the  structural  idealization  was  intended  to 
simulate  the  restraining  effect  of  the  epoxy  bond  between  the  waveguide  walls  and 
the  dielectric  material.  The  technique  by  which  the  spring  constant  k.,  was 
quantified  is  outlined  below. 

The  cross  sections  of  a  propellant  strand  and  a  casting  resin  section  were 

identical  and  recta  "ular.  However,  in  order  to  keep  the  input  to  the  finite  element 

model  computationally  simple  and  to  slant  the  analysis  toward  a  less  stiff  system 

(i.c..  a  more  conservative  system),  the  cross  sectional  areas  were  considered  circular 

and  were  constrained  to  have  the  same  area  as  the  rectangular  waveguide  cross 
sections.  Spring  constant  k,  was  determined  by  applying  circular  Hat  plate  theory 
relevant  to  deflections  small  in  comparison  to  plate  thickness. 
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If  a  material  obeys  Hooke’s  law,  the  shearing  stress  r  and  shearing  strain  7 
are  related  by(24) 


7  = 


T 

G 


(D-l) 


where  G  is  the  shear  modulus,  G  =  E/2(l+i>).  For  small  shearing  strains,  the 
following  approximation  can  be  written 


dw  r 

tan  7  =  —  =  7  =  —  (D-2) 

dr  G 


Define  a  spring  constant  given  by 


k  = 


(fs)r=R  o 
w 


(D-3) 


where  (fs)r=Ro  is  the  shear  force  at  the  edge  of  a  circular  flat  plate  of  thickness  h 
and  radius  Ro;  w  is  the  mean  deflection.  This  mean  deflection  can  be  expressed  as 


w 


/a  w  dA 
'Ro 


(D-4) 


By  the  use  of  equations  (D-l)  and  (D-2)  the  following  relation  can  be  written 


(f,)r- 


(D-5) 


l 
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The  spring  constant  of  equation  (D-3)  can  now  be  expressed  as 


(D-6) 


where  w  represents  the  local  deflection  of  the  circular  flat  plate  of  thickness  h. 
Assume  that  tire  plate  is  simply  supported  along  its  outer  edge  and  that  it  is 
subjected  to  a  load  per  unit  area,  q,  distributed  uniformly  over  its  area.  The 
expression* 25 *  for  the  local  deflection  of  such  a  plate  is 


w 


64  D 


(R:  - 


r2) 


/  5  +  v  ■. 

- R2 

\  1  +  v  ° 


n!ii  {*  +  v)  (R2 

8D  ( 1  -  v2 )  ° 


-r2) 


(D-7) 


where  v  is  Poisson’s  ratio  and  l)  is  the  flexural  rigidity  given  by 
D  =  Eh3 /1 2(1  -  v2 ).  The  second  term  on  the  right-hand-side  of  equation  (D-7) 
represents  the  correction  for  shearing  stresses  and  lateral  pressure. 

Substitution  of  equation  (D-7)  into  equation  (D-6)  and  subsequent 
evaluation  of  the  resulting  expression  yields  a  spring  constant.  Assume  that  this 
spring  constant  can  be  used  to  represent  the  spring  constant  k,  shown  in  Figure  3. 
To  compute  the  value  of  k-,  relevant  to  a  given  plate  material,  values  for  E,  h.  v. 
and  R<(  are  required.  The  elastic  modulus  tor  M6  propellant  is  114,000  psi.*26)  It 
was  assumed  that  a  Poisson’s  ratio  of  0.45  was  applicable  to  M6  propellant.  The 
cross  sectional  area  of  the  rectangular  waveguide  housing  the  dielectric  material  was 
0.1186  in.2;  hence,  the  radius  of  a  circular  cross  section  having  the  same  area  is 
0.1943  in.  Since  the  dielectric  material  was  divided  into  50  segments,  the  value  for 
h  was  0.10132  in.  Based  upon  the  above  parameter  values,  k,  was  computed  to  be 
3.6693  X  10s  lbf/in. 
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The  transient  response  analysis  for  the  structural  idealization  of  Figure  3 
was  achieved  via  use  of  a  finite  element  model.* 2 7)  The  time  dependent  loading  was 
considered  to  be  applied  to  the  first  interconnection  node.  In  the  model,  initial 
conditions  are  specified  with  respect  to  an  initial  equilibrium  configuration.  Attention 
was  focused  on  that  output  of  the  numerical  solution  which  yielded  the 

displacement,  stress,  and  velocity  of  the  first  interconnection  node.  Although  the 
model  allows  damping,  no  damping  was  considered  because  of  the  difficulty  in 

quantitatively  defining  the  inass  proportional  and/or  stiffness  proportional  damping 
constants  relevant  to  the  structural  idealization.  Hence,  only  undamped  transient 
numerical  solutions  were  analyzed. 
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I 

I 

i 
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PHASE  ANGLE  DERIVATIVE  RELATIONSHIPS 


After  noting  that  A,  a,  Bo,  and  (i  are  constants  while  t/>a,  ,  and  C  are 

functions  of  time,  take  the  time  derivative  of  equation  (16).  Multiply  the  result  by 
e'‘^a  and  then  equate  the  real  and  imaginary  parts  to  get 


•  jU  * 

p^/aBoe  a  +  C*m  sin(^m  -  4>a)  =  C  cos(\pm  -  \pa)  (E-l) 

•  u  6  • 

\j/  B  e  11  -  Ct p  .  cos  ( i//  -  4/  )  =  C  sin  ( i//  -  ^  )  (E-2) 

.1  <>  Mn  VTin  ra'  ’ni  ra  ' 


Dividing  equation  (E-2)  by  equation  (E-l)  and  then  multiplying  each  side  of  the 
resulting  equation  by  cos(\!/m  -  \!>a)  yields 


(E-3) 


If  after  squaring  equations  (13)  and  (14)  they  are  added  together,  the  resulting 
expression  becomes 


Be 


a  cos  ( ip m  -i J/a)  =  C  - 


A  cos  ( i p  -  a) 


(E-4) 


1 


209 


By  eliminating  C  from  equations  (13-15),  the  following  result  can  be  written: 


sin  (i//m  -  0a)  ~  -  A  sin  (0m  -  a) 


(E-5) 


Substituting  equations  (H-4)  and  (E-5)  into  equation  (E-3)  gives  the  desired  result: 


0,  = 


a  A  juA 

1-  -cos(0m  -  or)+  —  sin  (0m  -  a) 


(E-6) 
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In 

experiments 

A. 


EQUIPMENT  LIST 


the  following  five  lists  the  equipment  relevant  to  the  microwave 
is  identified. 

Electronic  components  and  signal  conditioning  devices  as  identified  in 
Figure  13: 

1.  Sweep  Oscillator  (microwave  source),  Alfred  Electronics,  Model 
650  main  frame  with  Model  655AK  oscillator  plug-in. 

2.  Synchronizer.  Frequency  Engineering  Laboratories,  Model  FEL 
133A. 

3.  RF  Power  Meter,  (monitor  of  microwave  source  output  power), 
Hewlett-Packard,  Model  43 IB. 

4.  Digital  Frequency  Meter  (10  Hz  resolution),  Hewlett-Packard, 
Model  5248M  main  frame  electronic  counter  with  Model  5255A 
plug-in  frequency  converter. 

5.  Microwave  Amplifier  (TWT),  Hewlett-Packard,  Model  495A. 

6.  Harmonic  Frequency  Converter,  Hewlett-Packard.  Model  84 11  A. 

7.  Network  Analyzer,  Hewlett-Packard,  Model  8410A. 

8.  Relative  Gain  Indicator.  Hewlett-Packard,  Model  8413A. 


9.  Amplifier  (20  dB),  Hewlett-Packard,  Model  465A. 
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10.  Filter  (passive  highpass  on  the  278  kHz  test  signal  entering  the 
phasemeter),  in-house  construction. 

11.  Phasemeter,  Dranetz  Engineering  Laboratories,  Inc.;  Model  305C, 
option  103,  Serial  No.  5010958405;  with  Plug-in  Model  305-PA-3002, 
Serial  No.  4121407005. 

12.  Filter-Amplifier  (16  lowpass  settings),  in-house  construction. 

13.  A/D  Data  Acquisition  System,  in  house  construction. 

a.  Memory:  three  channels  each  with  storage  capability  of  2048 
1 2-bit  words  and  one  channel  with  storage  capability  of 
2048  12  or  16-bit  words;  memory  construction  from  1024 
word  by  one  bit  static  random  access  memory  elements; 
these  Model  2102  RAM’s  were  manufactured  by  Fairchild 
Semiconductor. 

b.  12-bit  channels:  sample-and-hold  amplifiers,  Analog  Devices, 
model  SHA-5;  A/D  converter,  Analog  Devices,  Model 
ADC12QU. 

c.  15-bit  channel:  complete  A/D  conversion  subsystem,  Phoenix 
Data  Inc.,  Model  No.  DAS  7215-001-04-02-1,  Serial  No. 
A6884. 

14.  Piezoelectric  Quartz  Pressure  Transducers  (two  items),  Kristal 
Instrument  Corporation,  Model  6202  with  a  thermoprotective 
plate  impregnated  with  low  viscosity  oil. 


15.  Charge  Amplifiers  (two  items),  Kistler  Instrument  Corp., 
Model  503 
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16.  Strip  Chart  Recorder,  Teehni-rite  Electronics,  Inc.,  Model  TR-888  main 
frame  with  4  Model  TS C-810  signal  conditioning  plug-ins. 

17.  Time  Delay  Sequencer,  in-house  construction. 

18.  Capacitive  Discharge  Circuit,  in-house  construction. 

19.  Audio  Cassette  Interface,  in-house  construction. 

20.  Audio  Cassette  Tape  Recorder,  Superscope,  Model  CD-302A. 

B.  Microwave  circuitry  components  as  shown  in  Figure  13;  coaxial  components 
have  type  N  connectors  while  all  other  components  have  X-band  waveguide 
flanges; 

1 .  Ferrite  Isolators,  9  items,  Systron-Donner,  Model  DBG-480A. 

2.  3  dB  Directional  Couplers.  3  items.  Systron-Donner,  Model 

DBG-675-3  (2  items  used  in  the  reversed  mode). 

3.  Coaxial  Hybrid  Junction  (used  to  equally  divide  the  input  power  into 
two  transmission  lines),  1  item,  Narda  Microwave  Corporation.  Model 
3035. 

4.  10  dB  Directional  Couplers,  3  items.  Systron-Donner,  Model 

DBG -675-10. 

5.  20  dB  Coaxial  Directional  Coupler.  1  item,  Narda  Microwave 
Corporation.  Model  3045C-20. 

6.  3  Port  Circulator,  2  items.  Systron-Donner,  Model  DBG-490. 
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7.  Variable  Attenuators  (used  in  the  main  arms),  2  items, 
Systron-Donner,  Model  DBG-430. 

8.  Variable  Attenuators  (used  in  cancellation  loops),  2  items,  Arra 
Inc.,  Model  X520-40A. 

9.  Variable  Phase  Shifter  (used  in  main  test  arm),  1  item, 

Hewlett-Packard,  Model  X885A. 

10.  Variable  Phase  Shifters  (used  in  cancellation  loops),  2  items, 

Systron-Donner.  Model  DBG-915-1. 

11.  E-H  Tuners,  2  items,  Waveline  Inc..  Model  659. 

C.  Electronic  equipment  not  specifically  shown  in  Figure  13  but  still 
utilized  during  the  performance  of  a  test: 

1 .  Digital  Power  Meter  (used  to  check  input  power  to  the  Harmonic 
Frequency  Converter),  Pacific  Measurements,  Model  1018A. 

2.  RMS  Voltmeter  (used  for  noise  measurements  and  to  check 
initial  test  signal  voltage  amplitude  to  the  phasemeter), 
Hewlett-Packard,  Model  3400A. 

3.  Oscilloscope  (used  for  general  trouble-shooting  and  noise 
measurements),  Tektronix,  Model  545A  with  Model  1A7A 
high-gain  differential  amplifier  plug-in. 

4.  Precision  Voltage  Source,  in-house  construction. 

5.  Digital  Multimeter,  Hewlett-Packard,  Model  3465A. 


I. 


215 


6.  Digital  Voltmeter,  Hewlett-Packard,  Model  3460B. 

D.  Major  components  of  the  pressurizing  system  shown  in  Figure  20: 

1.  Gas  Booster  Pump  (up  to  40,000  psi),  American  Instrument 
Company,  Model  46-13475. 

2.  Pneumatically  Operated  Valves  (three  items),  American  Instrument 
Company,  Model  44-19731. 

3.  Check  Valves  (two  items),  American  Instrument  Company,  Model 
44-16380. 

4.  Dial  Pressure  Gages  (three  items:  used  for  monitoring  pressure  in 
booster  pump  and  in  combustion  fixture  chambers;  Bourdon 
type).  American  Instrument  Company.  Model  47-18350. 

5.  Pressure  Regulator.  Victor  Equipment  Company.  Model  SR-200C. 

E.  Combustion  fixture  sealing  components  relevant  to  Figures  18  and  19: 

1.  Self-Sealing  Electrical  Feedthroughs  (two  items),  American 
Instrument  Company.  Model  45-17582. 

2.  Rupture  Disc  Assembly  (thrust  nut  was  greatly  modified  before 
use),  American  Instrument  Company, 

a.  Thrust  nut.  Model  1636-147-1 


b.  Thrust  collar.  Model  1620-090 
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c.  Rupture  disks  (at  a  variety  of  failure  pressures),  Model 
45-19290-1 

3.  Major  End  Cap  Face  Seals 

a.  Resilient  metal  c-seals  (two  items)  Pressure  Science  Inc., 
Model  615A5X-0088-1. 

b.  Elastomeric  back-up  seals  (two  items),  Parker  Seal  Company, 
Model  2-(362)N674-70. 

4.  Flow  Nozzle  Face  Seal  (elastomeric),  Parker  Seal  Company, 
Model  2-(224)N674-70. 

5.  Explosive  Valve  Face  Seal  (resilient  metal  c-seal),  Pressure  Science 
Inc.,  Model  6 14A5 1-0050-1. 

6.  Dielectric  Transition  Window  Face  Seals 

a.  Larger  resilient  metal  c-seal.  Pressure  Science  Inc.,  Model 
614A5X-0044-1 . 

b.  Smaller  resilient  metal  c-seal.  Pressure  Science  Inc.,  Model 
633A5 1-0025-1. 
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DIGITAL  SIGNAL  ANALYSIS 


Digital  signal  analysis  was  performed  on  all  data  acquired  by  the  A/D 
system.  The  purpose  of  this  analysis  was  to  eliminate  unwanted  high  frequencies 
which  carried  little  or  no  information  and  to  produce  an  output  set  with  highly 

reduced  noise.  The  first  operation  done  on  a  given  discrete  sample  set  was  the 
computation  of  the  spectrum  which  yielded  the  distributions  over  frequency  of  both 
the  amplitude  and  the  energy.  To  calculate  these  distributions,  the  fast  Fourier 
transform  (FFT)  was  utilized.  This  transformation  is  a  very  efficient  algorithm  for 
computing  the  finite  discrete  Fourier  transform  (DFT).<38'40>  For  the  FFT 
algorithm  used  with  the  experimental  data,  the  number  of  samples  which  could  be 
operated  upon  was  constrained  to  be  an  integer  power  of  2. 

Of  course,  the  experimental  data  of  interest  did  not  always  have  precisely 
2m  samples  in  a  set.  However,  the  actual  number  of  data  samples  in  a  time  interval 
of  interest  was  always  just  slightly  more  or  somewhat  less  than  a  particular  2m 
value.  Hence,  for  each  physical  quantity  being  analyzed,  the  most  appropriate  2m 

number  of  discrete  values  was  operated  upon  by  the  FFT  algorithm.  The  resulting 
distortion  in  the  spectral  information  was  slight  and  did  not  compromise  the 

judgment  of  where  an  acceptable  digital  lowpass  cutoff  frequency  could  be 

established  for  a  given  transient  signal.  It  should  be  noted  that  all  spectral  analysis 
was  done  upon  physically  dimensioned  discrete  data  which  had  neen  amplitude 
translated  such  that  amplitude  changes  began  from  zero  values. 

The  Fourier  transform  of  any  real  function  fit)  is 

F(joo)  =  f  fitje^'dt  (G-l ) 
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The  DPT.  denoted  Fju>).  can  be  defined  using  the  zero-order  approximation  to 
F(jcj):  that  is. 


F(joa)  =  £  f(nT)e  iwnT 

n  -  - 00 


(G-2) 


where  T  is  the  sampling  interval  and  also  the  interval  of  integration  in  the  time 
domain.  An  alternate  way  of  describing  the  evolution  of  the  DFT  is  to  note  that  the 
zero  order  approximation  of  F(jto)  divided  by  T  equals  F(joo).  Consequently, 

✓ 


TF(jco)  =  F(jco) 


(G-3) 


Any  actual  computation  of  the  DFT  must  involve  a  finite  summation  of  terms; 
hence,  suppose  there  are  N  samples  of  f(t).  The  DFT  is  periodic  over  cj  with 
period  27r/T.  Hence,  by  letting 

2  rnn 

w  =  -  ;  in  =  0.1.  ...  .  N  -  1  (G-4) 

NT 


the  usual  formula  lor  computing  the  independent  values  of  the  DFT  becomes 


N  I  j(2»mn/N) 

1  !ne  ;  m 

n  l» 


0.1.  ....  N  -  1 


(G-5) 


» 
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where  f(nT)  has  been  written  as  f'n.  Actually,  only  the  values  of  Fm  for  m  =  0  to 
N/2  really  need  to  be  computed. 

For  each  transient  data  signal,  approximations  were  determined  for  |F(jtu)| 
and  |F(jcj)|2,  the  Fourier  amplitude  and  energy  spectrums,  respectively.  These 
approximations  were  T(F(jco)l  and  T2  |F(jcu)l2 ;  they  required  the  computation  of  the 
DFT  of  a  given  sample  set  [  f n  ]  acquired  from  the  A/D  system.  It  was,  of  course, 
assumed  that  the  sampling  interval  T  had  been  chosen  so  that  essentially  all  of  the 
spectral  content  of  the  waveform  was  contained  below  1/2T  Hz,  the  Nyquist  or 
folding  frequency.  The  DFT  amplitude  was  found  from  the  quadrature  of  the  real 
and  imaginary  parts;  i.e.. 


IF„ 


/ 

N  - 1  27rmn 

2 

N-i  2rrmn 

I 

2  f  cos 

+ 

2  f  sin 

V 

n  =  0  "  N 

n  =  0  "  N 

2y/2 


(G-6) 


where  m  needed  only  to  vary  from  0  to  N/2. 

Instead  of  utilizing  T|F(jw)l  or  T2  |F(jw)l2  directly,  two  hybrid  expressions 
were  constructed  and  plotted  as  a  function  of  frequency;  the  frequency  range  was 
v  =  0  to  1/2T  and  the  increment  was  =  1/NT.  The  expressions  for  relative 
amplitude  and  accumulative  energy  are 

(G-7) 

(G-8) 


i 
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In  equation  (G-7),  the  term  |F  |  is  simply  the  maximum  value  of  the  DFT 
amplitude.  Hence,  the  plot  of  Am  versus  v  showed  the  spectrum  in  decibels  below 
the  maximum  |Fm  |.  Values  of  E  can  range  from  0  to  1  and  they  have 

significance  only  if  essentially  all  the  spectral  content  of  the  transient  signal  is 
contained  below  1/2T  Hz.  By  analysis  of  the  spectral  information  contained  in  the 
plots  of  A  and  E  versus  v.  an  appropriate  digital  lowpass  cutoff  frequency  was 
established  for  each  transient  signal  acquired  by  the  A/D  system.  A  cutoff  frequency 
was  chosen  where  A)n  was  sufficiently  low  and  E  sufficiently  high.  Generally, 
these  threshold  values  were  A  =  -  40  dB  or  less  and  E  =  0.995  or  more. 

m  m 

Once  a  decision  had  been  made  on  what  frequencies  in  the  signal  to  pass 
and  which  to  reject,  the  given  sample  set  [f  |  was  passed  through  a  lowpass, 

inphase  nonrecursive  digital  filter  which  then  created  the  output  set  [gm].(4l)  A 
nonrecursive  filter  is  one  in  which  the  formula  for  gm  contains  explicitly  only 
values  from  the  set  (fml.  A  general  form  of  the  linear  nonrecursive  algorithm  is 

Sn,  =  2  Nbnfm-n  (C-9) 

n  =  -  N 


where  the  limit  N  is  finite  to  insure  realization.  Note  that  there  are  a  total  of 
2N  +  1  filter  weights.  Each  coefficient  b  can  be  any  real  value  including  zero.  The 
specific  values  of  the  bn’s  are  dependent  upon  the  desired  real  or  complex  transfer 
function  defined  in  the  frequency  space.  The  counterpart  or  analog  to  the  set  of 
discrete  smoothing  weights  |bnl  is  the  continuous  weighting  function  b(t)  which 
would  be  associated  with  continuous  input  and  output  functions;  that  is,  by 
convolution 


g(t)  = 


f(t  -  T)b(T)dr  = 


f(r)b(t  -  T)dr 


(G-10) 
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The  bn  can  be  evaluated  directly  from  b(t)  by  noting  that  bn  =s  b(nT)T  where 
equidistant  samples  have  been  assumed. 

Assume  that  the  transfer  function  H(co)  is  real  and  that  it  is  symmetric 
with  respect  to  co  =  0;  i.e.,  H(co)  =  H(-co).  The  weighting  function  b(t)  is 

c  00 

b(t)  =  -  \  H(cc)ejuJ  ‘dco  =  — 

2ir  1  2n 


£ 


‘dw 


(G-l  1 ) 


Til  is  equation  can  be  written  in  real  form  over  the  positive  frequency  space.  The 
result  is 


b(t)  = 


_l_ 

it 


H(co)  cos  cot  dco 


(G-l  2) 


where  b(t)  is  real  and  is  symmetric  with  respect  to  t  =  0;  i.e..  b(t)  =  b(-t).  If  the 
time  derivative  of  the  output  function  is  desired,  then  from  equation  (G-10) 


g(t)  = 


d[g(t)| 

at 


f( t )  —  [b(t  -  r)jdr 

at 


f(r)c(t  -  r)dr 


(G-13) 


where  c(t')  =  a/at'  [b(t')|  when  t'  =  t  -  t.  By  the  use  of  equation  (G- 11).  c(t)  can 
be  expressed  as 


c(t) 


2n 


jcoH(oj)eiw  ’dco 


!-  jcoH(co)e  ^'dco 

00 


(G-14) 
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Clearly,  jcoH(co)  is  the  Fourier  transform  of  c(t)  while  H(co)  is  the  Fourier 
transform  of  b(t).  Rewriting  equation  (G-14)  in  real  form  over  the  positive 
frequency  space  yields 


c(t) 


coH(co)  sin  cot  dco 


CG- 15) 


where  c(t)  is  a  real  function  and  is  antisymmetric  with  respect  to  t  =  0;  i.e., 
c(-t)  =  -c(t). 

The  specific  form  of  the  transfer  function  H(co)  used  in  the  digital  signal 
analysis  of  the  experimental  data  is  shown  in  Figure  45.  This  transfer  function  gave 
a  filter  having  unity  gain  and  zero  phase  shift  in  the  pass  band  with  a  first  order 
roll-off  between  the  pass  (go  <  eot)  and  rejection  (go  >  co2)  bands.  After  the  direct 
application  of  equations  (G- 12)  and  (G-15),  the  smoothing  and  derivative  weighting 
functions  for  this  particular  H(co)  were  computed  to  be 


b(t)  = 


cos  co  |  t  -  cos  co2t 
tt(co2  -  co,)t2 


(G-16) 


c(t) 


co2 1  sin  co2 1  -  coj  t  sin  co,  t  +  2  cos  to2 1  -  2  cos  tO]  t 
7T(C02  -  CO,  )t3 


(G-l  7) 


The  discrete  weights  needed  in  the  algorithm  represented  by  equation  (G-9),  come 
directly  from  equations  (G-16)  and  (G-l 7).  These  discrete  weights  are 


b 


S 

n 


Tb(nT)  = 


cos(27rf,  nT)  -  cos(2irf2nT) 
27r2n2  T  (f2  -  f,) 


(G-l  8) 
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b 


d 

n 


Tc(nT)  = 


f2sin(27rf2nT)  -  ft  sin(2ffflnT) 
7rTn2  (f2  -  f, ) 


(G-19) 


The  superscripts  “s”  and  “d”  refer,  respectively,  to  smoothing  and  derivative. 
Because  of  the  symmetry  of  b*  and  the  antisymmetry  of  bj|,  it  was  unnecessary  to 
compute  these  coefficients  for  negative  n  values.  The  specific  algorithm  used  to  get 
the  smoothed  output  was 


a  = 

em 


N 

2 

n  =  -  N 


bsf 

n  ni 


=  bs 

n 


b0  = 


T(f2  +  f, ) 


(G-20) 


To  compute  the  derivatives  of  the  output,  the  following  algorithm  was  utilized 


g 


d 

m 


b  f 

n  m  -  n 


-  bj  ;  1,;  •  0 


(G-2 1 ) 


The  lowpass  transfer  function  shown  in  Figure  45  is  a  Fourier  transform 
associated  directly  with  b(t)  and  indirectly  with  c(t).  For  n  in  the  interval  [-N,N], 
the  discrete  weights  b*  and  bjj  are  associated  with  a  finite  DFT  which  is  a 
least-squares  approximation  to  H(co).  With  N  finite,  the  goodness  of  fit  between 
H(oj)  and  H(co)  depends  upon  the  sharpness  of  the  roll-off  after  co,  together  with 
N,  the  parameter  which  defines  the  number  of  weights.  Central  to  this  topic  of 
error  between  the  desired  and  the  realized  spectral  window  is  the  so-called  Gibbs 
phenomenon  which  manifests  itself  as  a  fixed  percentage  overshoot  and  ripple  before 
and  after  an  approximated  discontinuity/4  2  * 


For  a  lowpass  digital  filter  having  a  desired  transfer  function  identical  to 
that  in  Figure  45,  Ormsby*4  3)  discussed  the  mean  square  error  e  between  H(oo)  and 
H(w).  His  empirical  results  for  e  <  0.006  and  e  <  0.015  were  correlated  with  the 
following  equation 


0.012 

XrN 


(G-22) 


where 


X' 


"  kXR 


(G-23) 


The  sharpness  of  the  rod-off  after  ocj  is  specified  by  XR  and  k  is  a  constant 
greater  than  one  which  defines  the  size  of  the  transition  band  in  which  the  effective 
error  is  e.  For  example  with  k  =  1.2,  X  =  co/gjs,  and  Xc  =  co,/cos,  the  maximum 
errors  are  less  than  e  for  X  <  (Xc  -  0.1  XR)  and  X  >  (X  +1.1  XR).  If  in 
equation  (G-22)  the  value  of  N  is  made  increasingly  large,  then  a  smaller  XR  can 
exist  for  the  desired  maximum  error  or  a  smaller  e  can  be  gained  for  the  same  XR. 

The  parameters  which  indicate  the  high  accuracy  or  low  error  in  the 
digital  filtering  conducted  on  the  data  from  the  steady-state  burning  tests  are  shown 
in  Table  5.  After  the  analysis  of  the  spectral  content  of  all  the  discrete  physical 
data  relevant  to  a  given  test,  a  single  digital  lowpass  cutoff  frequency  was  found  to 
suffice  for  all  transient  signals.  Each  value  of  f,  in  Table  5  represents  such  a 
frequency.  In  order  to  justify  the  choice  of  N  =  200  as  shown  in  Table  5.  consider 
the  following  argument.  Let  e  =  0.006  and  N  =  49,  then  from  equation  (G-22)  XR 
is  found  to  be  0.0408.  Since  XR  is  0.04,  k  will  be  1.02.  Clearly,  by  increasing  N 
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from  49  to  200  the  error  e  will  get  smaller  for  the  same  XR  ■  Hence,  by  choosing 
N  =  200  the  smoothed  output  will  have  essentially  no  error.  Also,  by  making  N 
larger  the  derivative  of  the  output  becomes  a  better  approximation  to  the 
desired  though  not  realizable  case  of  N  =  <*>.  It  is  important  to  note  that  for  equal 
X'R  the  error  in  is  much  more  sensitive  to  N  than  is  the  error  in  g^ .  Thus, 
when  derivatives  are  to  be  calculated  with  the  differentiating  digital  filter  represented 
by  equations  (G-19)  and  (G-21),  the  value  chosen  for  N  should  be  larger  than  the 
number  which  would  be  satisfactory  if  only  g^  were  being  computed.  This  is 
another  reason  why  N  was  made  200  instead  of  only  49.  For  all  steady-state 
burning  tests,  time  derivatives  were  calculated  by  the  digital  nonrecursive 
differentiating  filter. 

Presented  in  Table  6  are  the  digital  filter  parameters  relevant  to  all 
dynamic  tests.  Unlike  the  steady-state  tests,  no  single  lowpass  cutoff  frequency  f, 
was  applied  to  all  physical  quantities  of  interest  on  a  given  dynamic  test.  Hence,  no 
values  f,  are  indicated  in  Table  6.  For  the  dynamic  tests,  these  cutoff  frequencies 
and  the  associated  physical  quantities  are  given  in  Table  15  which  is  discussed  in 
Section  6.31  and  not  here.  As  indicated  by  the  values  of  XR  in  Table  6,  the 
roll-off  after  a  cutoff  frequency  was  sharper  than  that  for  a  steady-state  test. 
Consider  the  consequences  of  this  fact.  If  k  =  1.02  as  in  the  steady-state  example, 
then  XR  =  0.0102  when  XR  =  0.010.  When  e  =  0.006  and  X'R  =  0.0102  the  value 
of  N  from  equation  (G-22)  becomes  196.  Obviously  the  use  of  N  =  200  gave  very 
good  results  for  the  dynamic  g^  although  the  filtered  output  was  not  quite  as  error 
free  as  that  for  the  steady-state  tests.  A  problem  arose  when  the  dynamic  g^  was 
computed  with  N  =  200.  It  was  found  that  to  get  a  high  quality  magnitude 
response  from  the  differentiating  nonrecursive  algorithm,  it  was  necessary  to  use 
values  of  N  greater  than  200.  Such  increases  in  the  number  of  filter  weights  caused 
the  computer  costs  to  become  prohibitive.  One  remedy  to  this  problem  would  have 
been  to  use  a  differentiating  digital  filter  with  weighting  coefficients  different  from 
those  given  by  equation  (G-19).  There  are  many  other  differentiating  filters  and  Kuo 
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and  Raise/4 4  *  provide  a  comparison  of  the  amplitude  responses  of  several  different 
designs.  The  solution  which  was  implemented  was  to  apply  central  differences  to  the 
already  digitally  filtered  data.  Thus,  for  all  dynamic  tests,  time  derivatives  of  any 
sample  set  were  computed  in  this  manner. 
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